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Special  contract  SPC-97-4017 

Planetary  Waves  and  Mesoscale  Disturbances 
in  the  Middle  and  Upper  Atmosphere. 

Final  report 

I.  Measurements  and  analysis. 

1.1.  Emission  rate  and  temperature  observations  by  MORTI. 

Mesopause  Oxygen  Rotational  Temperature  Imager  (MORTI)  is  ground-based  optical  instru¬ 
ment  for  emission  rate/temperature  measurements  by  using  the  molecular  oxygen  nightglow 
emission.  Its  construction  was  described  in  details  by  Wiens  et  al.  (1991).  Here  we  present 
a  general  concept  of  an  instrument  construction  and  discuss  some  preliminary  results  of  observa¬ 
tions  in  Almaty.  MORTI  uses  the  (0-1)  band  of  O2  spectrum  around  866.5  nm.  The  instrument 
consists  of  a  conical  mirror,  a  Fresnel  lens,  an  interference  filter,  a  photographic  lens,  and  CCD 
camera.  The  optical  axis  of  the  MORTI  is  directed  vertically  and  the  conical  mirrpr  provides 
a  view  of  the  annulus  with  an  inner  radius  of  47  km  and  an  outer  radius  of  63  km  at  altitude  of 
95  km.  The  interference  filter  is  centered  at  867.6  nm  and  has  bandpass  of  0.27  nm.  It  transmits 
the  O2  lines  with  different  wavelength  at  different  angular  distances  from  the  optical  axis  form¬ 
ing  six  concentric  rings  at  the  CCD  plane  corresponding  to  six  pairs  of  the  rotational  lines.  CCD 
(charge  coupled  device)  camera  operates  as  a  detector  of  the  sky  image.  The  radial  distribution 
of  the  image  intensity  contains  the  information  on  a  rotational  temperature,  and  integration  of 
the  intensity  radially  over  the  image  gives  the  emission  rate  for  the  corresponding  sky  region. 
Evaluations  show  that  the  emission  rate/temperature  measurements  have  a  precision  of  ±  2% 
and  ±  2  K,  respectively.  Passing  overhead  acoustic-gravity  waves  (AGW)  modulate  the  rota¬ 
tional  temperature  and  the  emission  rate  so  the  temporal  behaviour  of  these  parameters  gives 
the  information  on  AGW  frequency  spectra  and  their  spatial  characteristics.  In  order  to  esti¬ 
mate  the  AGW  frequency  spectra  the  temperature  and  emission  rate  are  processed  by  a  Fourier 
analysis.  To  obtain  information  on  an  azimuth  and  a  phase  speed  of  the  AGW  propagation 
the  sky  field  of  view  is  divided  into  12  sectors,  and  time  series  of  the  temperature  and  emission 
rate  in  these  sectors  have  to  be  used  for  estimation  of  the  AGW  propagation  characteristics. 
The  apparent  phase  differences  between  similar  signatures  from  the  various  image  sectors  give 
the  information  on  the  azimuth  of  AGW  propagation  and  phase  speed. 

The  interference  spectro-photometer  MORTI,  using  advanced  CCD  imaging  technology,  has 
been  installed  at  the  Institute  of  Ionosphere  in  September  1997.  Since  October  1997  until  April 
1998  the  compaign  type  measurements  (10-12  days  a  month)  were  carried  out  during  clear 
sky  nights  when  the  moon  was  not  full.  To-  obtain  the  images  we  used  5  min  exposure  time. 
A  distance  of  10  km  between  Almaty  city  and  Institute  of  Ionosphere  ensure  good  conditions 
and  effectively  abolished  the  background  illumination. 


The  following  tables  show  the  statistics  of  our  observations 


October  1997 

November  1997 

23  (321  min) 

clouds 

1  (695  min) 

clear  sky 

24  (376  min) 

clouds 

2  (697  min) 

clouds 

25  (432  min) 

clouds 

3  (700  min) 

clear  sky 

26  (681  min) 

clear  sky 

4  (702  min) 

clouds 

27  (683  min) 

clear  sky 

5  (704  min) 

clouds 

28  (685  min) 

clear  sky 

6  (520  min) 

clouds 

29  (688  min) 

clear  sky 

7  (461  min) 

clouds 

30  (690  min) 

clear  sky 

8  (398  min) 

clouds 

31  (693  min) 

clear  sky 

9  (333  min) 

clouds 

November  1997 

December  1997 

21  (331  min) 

clouds 

1  (751  min) 

clouds 

22  (387  min) . 

clouds 

2  (752  min) 

clouds 

23  (442  min) 

clouds 

3  (753  min) 

clouds 

24  (498  min) 

clear  sky 

4  (754  min) 

clouds 

25  (743  min) 

clear  sky 

5  (556  min) 

clear  sky 

26  (744  min) 

clear  sky 

6  (492  min) 

clear  sky 

27  (745  min) 

clear  sky 

7  (427  min) 

clear  sky 

28  (747  min) 

clear  sky 

8  (361  min) 

clear  sky 

29  (749  min) 

clear  sky 

30  (750  min) 

clear  sky 

December  1997 

January  1998 

21  (359  min) 

clouds 

1  (761  min) 

clouds 

22  (413  min) 

clouds 

2  (760  min) 

clear  sky 

23  (468  min) 

clouds 

3  (759  min) 

clouds 

24  (523  min) 

clear  sky 

4  (511  min) 

clouds 

25  (764  min) 

clouds 

5  (444  min) 

clouds 

26  (763  min) 

clouds 

6  (378  min) 

clouds 

27  (763  min) 

clouds 

28  (763  min) 

clouds 

29  (762  min) 

clouds 

30  (762  min) 

clouds 

31  (762  min) 

clear  sky 

January  1998 

February  1998 

21  (413  min) 

clear  sky 

1  (716  min) 

clouds 

22  (469  min) 

clouds 

2  (496  min) 

clear  sky 

23  (525  min) 

clear  sky 

3  (429  min) 

clear  sky 

24  (581  min) 

clouds 

4  (361  min) 

clear  sky 

25  (730  min) 

clear  sky 

26  (728  min) 

clouds 

27  (726  min) 

clouds 

28  (724  min) 

clouds 

29  (722  min) 

clear  sky 

30  (720  min) 

clear  sky 

31  (718  min) 

clear  sky 

February  1998 

March  1998 

19  (362  min) 

clear  sky 

1  (646  min) 

clouds 

20  (418  min) 

clouds 

2  (643  min) 

clouds 

21  (473  min) 

clouds 

3  (640  min) 

clouds 

22  (526  min) 

clouds 

4  (383  min) 

clouds 

23  (663  min) 

clouds 

24  (660  min) 

clear  sky 

25  (657  min) 

clouds 

26  (654  min) 

clear  sky 

27  (652  min) 

clear  sky 

28  (649  min) 

clear  sky 

March  1998 

April  1998 

19  (266  min) 

clear  sky 

1  (552  min) 

clouds 

20  (321  min) 

clouds 

2  (290  min) 

clear  sky 

21  (373  min) 

clouds 

3  (269  min) 

clear  sky 

22  (423  min) 

clouds 

23  (471  min) 

clear  sky 

24  (577  min) 

clouds 

25  (574  min) 

clouds 

26  (571  min) 

clouds 

27  (568  min) 

clouds 

28  (565  min) 

clouds 

29  (562  min) 

clouds 

- 

30  (558  min) 

clear  sky 

31  (555  min) 

clouds 

4 


1.2.  Periodogram  analysis  of  MORTI  data. 

Lomb-Scargle  (LS)  periodogram  analysis  (Lomb,  1976;  Scargle,  1982)  was  first  applied  to  as¬ 
sess  the  statistically  significant  oscillations  present  in  the  data  (Press  and  Teukolsky,  1988). 
A  quadratic  fit  was  subtracted  from  each  data  set  to  remove  trends  whose  periods  were  com¬ 
mensurate  with  the  data  set  length.  This  was  necessary  to  minimize  spectral  leakage  and 
to  preserve  the  probability  distribution  of  the  periodogram  (Scargle,  1982;  Horne  and  Bali- 
unas,  1986).  The  typical  examples  of  the  data  observed  over  Almaty  (43.6  N,  76.6  E)  and 
the  results  of  the  periodogram  analysis  are  shown  in  Figs.  1-8.  Two  groups  of  statistically 
significant  (>  99%  confidence  level)  peaks  appear  in  Figs.  1-8.  The  lower  frequency  group  cor¬ 
responds  to  the  periods  of  5-9  hours  (Figs.  1-4,  7,  8).  These  variations  are  most  likely  due  to 
the  tidal  (6h  and/or  8h)  oscillations  similar  to  reported  by  Wiens  et  al.  (1995)  and  Ozonovich 
et  af.  (1997).  The  higer  frequency  group  corresponds  to  the  periods  of  2-3  hours  (Figs.  5,  6), 
and  these  variations  are  due  to  the  AGW  propagation  at  the  mesopause  heghts. 


Krassovsky’s  ratio  r)  (1972)  is  a  complex  quantity  which  can  be  considered  as  a  transfer  function 
between  planetary  wave,  tidal,  or  AGW  driven  temperature  oscillations  and  resultant  emission 
rate  oscillations.  In  a  ground-based  measurements  the  parameter  7]  may  be  expressed  as  follows 
(Zhang  et  ah,  1993b) 

OO 


B'/Bo 

TJTo 


f  e'dz 
0 _ 

OO  5 

feoV/Todz 

0 


(1) 


where  B  is  the  apparent  emission  rate,  Tm  is  the  measured  temperature  weighted  by  the  appar¬ 
ent  emission  rate,  and  e  is  the  volume  emission  rate.  The  subscript  zero  denotes  the  unperturbed 
values  and  the  perturbations  are  denoted  by  a  prime.  The  phase  of  Krassovsky’s  ratio  is  given 
by  relative  phase  of  the  emission  rate  and  temperature  oscillations  ^pri  =  ‘fB  ~  ^t-  Confusion 
exists  in  regard  to  the  sign  of  the  phase  of  Krassovsky’s  ratio  For  instance,  in  the  paper  Zhang 
et  al.  (1993a)  on  can  read  that  ’’when  >  0,  the  temperature  leads  the  emission  rate;  when 
(^,,  <  0  the  eihission  rate  leads”  (see  page  369),  however,  when  the  authors  discuss  the  phase 
of  Krassovsky’s  ratio  predicted  by  the  theory  (for  most  cases  the  phase  is  negative),  they  make 
a  conclusion  that  ’’the  temperature  fluctuations  lead  the  column  emission  rate  fluctuations  for 
all  wavelengths  and  periods”  (Zhang  et  al.,  1993b,  page  388).  It  seems  that  source  of  this 
ambiguity  is  a  possibility  to  determine  the  phase  of  Fourier  components  in  time  series  by  two 
different  ways:  as  the  cosine  argument  of  y{t)  at  the  time  i=0,  where  y{t)  =  Acos{u:t  -f  ip), 
A  is  amplitude,  and  tp  is  the  phase  (Hocke,  1998);  or  as  a  position  of  the  maximum,  i.e., 
y{t)  =  Acos{u}t  —  (p)  (Ozonovich  et  al.,  1997).  Another  possible  source  of  the  confusion  is 
the  different  choice  of  sign  in  the  presentation  of  perturbations  by  plane  waves.  To  avoid  the 
ambiguity  connected  with  the  sign  of  the  phase  of  Krassovsky’s  ratio  we  adopt  the  following 


26  November,  1997 


Universal  Time  (hr) 


5 


Figure  1:  The  relative  emission  rate  (a)  and  rotational  temperature  (b).  The  normalized  LS 
periodograms  (c)  from  S'/^o  -  solid  line  and  TI^/Tq  -  dashed  line  for  the  data  observed  over 
Almaty  on  Nov.  26  to  Nov.  27,  1997.  The  quadratic  trends  are  removed. 


ssion  ra 


Figure  2:  As  in  Fig.  1,  but  for  the  data  observed  on  Nov.  27  to  Nov.  28,  1997. 


Figure  3:  As  in  Fig.  1,  but  for  the  data  observed  on  Nov.  28  to  Nov.  29,  1997. 


2  January,  1998 


8 


Period  (min) 


Figure  4:  As  in  Fig.  1,  but  for  the  data  observed  on  Jan.  2  to  Jan.  3,  1998 


Lmi 


2  February,  1998 


Figure  6:  As  in  Fig.  1,  but  for  the  data  observed  on  Feb.  2  to  Feb.  3,  1998 


Temperature 


26  February,  1998 


Figure  7;  As  in  Fig.  1,  but  for  the  data  observed  on  Feb.  26  to  Feb.  27,  1998 
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presentation  of  perturbed  quantities 


(2) 


where  the  phase  <^(2)  is  the  argument  of  complex  amplitude  fkx,uj{z)-  It  is  possible  to  show 
that  such  determination  of  the  phase  of  the  wave  is  equivalent  to  the  second  determination  of 
the  phase  in  Fourier  analysis  of  time  series,  i.e.,  y{t)  =  Acos{u)t  —  (p).  In  our  determination  the 
phase  of  Krassovsky’s  ratio  is  positive  if  the  temperature  leads  the  emission  rate,  and  negative 
if  the  emission  rate  leads  the  temperature.  It  is  necessary  also  to  note  that  our  agreement  is 
different  from  the  one  traditionally  accepted  in  the  papers  on  airglow  emissions,  but  usually  it 
is  accepted  in  the  hydrodynamics. 


In  order  to  obtain  the  amplitude  and  phase  of  Krassovsky’s  ratio  from  the  experimental  data, 
we  used  the  LS  periodogram  method  for  estimation  the  amplitude  and  phase  spectra  (Hocke, 
1998).  The  results  obtained  with  data  shown  in  Fig.  1-8  are  presented  in  Table  1. 


Table  1 


Date 

Period  (min) 

B'fBo  (percent) 

T'/Tq  (percent) 

1^1 

(deg) 

26.11.97 

377.3 

23.7 

4.2 

5.7 

37.8 

198.0 

5.9 

1.6 

3.6 

30.1 

102.5 

2.6 

Uncert. 

Uncert. 

Uncert 

27.11.97 

541.5 

20.2 

4.8 

4.2 

16.6 

297.8 

8.2 

Uncert. 

Uncert. 

Uncert. 

165.4 

6.5 

1.5 

4.4 

17.9 

28.11.97 

406.4 

27.8 

3.9 

7.1 

21.8 

183.5 

5.0 

Uncert. 

Uncert. 

Uncert. 

127.9 

3.3 

Uncert. 

Uncert. 

Uncert. 

2.01.98 

479.2 

18.4 

3.0 

6.1 

10.3 

250.0 

3.8  • 

Uncert. 

Uncert. 

Uncert. 

166.5 

3.7 

1.15 

3.2 

2.7 

21.01.98 

142.8 

6.3 

1.9 

3.3 

32.2 

100.3 

2.7 

Uncert. 

Uncert. 

Uncert. 

2.02.98 

149.8 

5.9 

1.6 

3.8 

22.1 

102.4 

3.5 

0.65 

5.3 

6.4 

26.02.98 

388.8 

35.5 

4.9 

7.3 

20.2 

206.4 

9.1 

1.9 

4.9 

11.3 

127.6 

4.4 

Uncert. 

Uncert. 

Uncert. 

83.6 

4.3 

Uncert. 

Uncert. 

Uncert. 

23.03.98 

295.5 

19.1 

5.8 

3.3 

19.7 

154.1 

8.1 

Uncert. 

Uncert. 

Uncert. 

99.7 

4.9 

Uncert. 

Uncert. 

Uncert. 

65.2 

4.3 

jjncert. 

Uncert. 

Uncert. 

54.7 

3.9 

Uncert. 

Uncert. 

Uncert. 
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This  table  shows  that  in  general  the  amplitude  of  the  emission  rate  and  temperature  oscilla¬ 
tions  decreases  with  the  decreasing  of  the  wave  period.  The  comparison  of  \t]\  and  (/?,,  obtained 
from  our  measurements  with  the  ones  obtained  in  other  investigations  (Zhang  et  ah,  1993a, b; 
Ozonovich  et  ah,  1997)  shows  a  good  agreement.  The  largest  amplitudes  of  Krassovsky’s  ra¬ 
tio  are  obtained  for  long-period  waves  (28.11.97  and  26.02.98),  but  clear  dependences  of  the 
amplitude  and  phase  of  t]  on  the  wave  period  are  not  evident.  Perhaps  the  later  is  due  to  dif¬ 
ferent  nature  of  long  period  waves  [migrating  and  nonmigrating  tides  or  inertia-gravity  waves, 
see  the  discussion  in  Ozonovich  et  al.  (1997)  and  Lindzen  (1984)].  In  all  cases  the  phase  of 
Krassovsky’s  ratio  is  positive  and  the  temperature  variations  leads  the  emission  rate  variations. 
Our  results  also  show  that  the  determination  of  the  wave  parameters  from  the  measurements 
of  the  emission  rate  variations  is  more  reliable  then  from  the  measurement  of  the  temperature, 
because  the  temperature  fluctuations  are  not  as  large  as  the  emission  rate  fluctuations  (the 
amplitude  of  Krassovsky’s  ratio  changes  from  3  to  7).  However,  it  is  necessary  to  note  that 
thin  tropospheric  clouds  can  change  the  measured  emission  rate  without  appreciably  changing 
the  measured  rotational  temperature,  because  we  see  no  reason  that  the  clouds  would  change 
the  spectral  distribution  from  which  we  determine  the  temperature  (Wiens  et  ah,  1991). 

Numerical  calculations  of  Krassovsky’s  ratio  were  performed  using  the  approach  considered 
in  Zhang  et  al.  (1993b).  The  AGW  propagation  at  the  mesopause  heights  was  simulated  with 
the  numerical  model  developed  by  Bidlingmayer  and  Pogoreltsev  (1992).  The  AGW  simulations 
were  performed  for  nonisothermal  atmosphere,  but  without  account  of  the  background  wind.  It 
is  possible  to  include  the  background  wind  into  consideration  (Pogoreltsev  and  Pertsev,  1995), 
but  we  have  not  yet  the  software  for  devision  of  images  into  12  sectors  and  we  cannot  determine 
the  azimuth  of  AGW  propagation  (Dr.  Rudy  Wiens  promissed  to  provide  us  with  this  software 
shortly).  The  calculations  were  carried  out  for  the  model  atmosphere  (Hedin  et  ah,  1996). 


Table  2  lists  theoretical  calculations  of  the  amplitude  and  phase  of  Krassovsky’s  ratio. 


Table  2 


Frequency 

rad/sec 

km 

1^1 

(deg) 

7  •  10-^ 
(T=150  min) 

600 

6.96 

54.8 

1000 

4.54 

38.3 

1400 

4.16 

27.4 

2.5  •  10-^ 
(T=400  min) 

2000 

5.87 

50.7 

3000 

4.53 

35.8 

4000 

4.20 

26.2 

5000 

4.17 

17.2 

One  can  see  that  the  results  of  simulation  are  reasonable  and  can  explain  the  most  of  ex¬ 
perimental  results  listed  in  Table  1. 
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1.3.  Cophase  analysis  of  MORTI  data. 

Although  we  cannot  apply  directly  the  sector  analysis  to  the  MORTI  data  obtained  in  Al¬ 
maty,  it  is  possible  to  develop  the  software  which  will  be  used  in  future,  when  we  will  have  the 
program  of  division  the  MORTI  image  into  12  sectors.  For  this  purpose  MORTI  data  obtained 
at  Bear  Lake  Observatory  in  January  1993  (Wiens  et  al.,  1997)  have  been  used.  The  MORTI 
image  maps  an  annular  field  of  view  down  to  a  solid  circle  in  which  azimuth  is  preserved.  By 
analyzing  each  image  in  12  sectors,  12  time  series  become  available  in  relative  temperature  and 
relative  column  emission  rate  for  fixed  point  of  sky.  The  amplitude  and  phase  spectra  can  be 
found  by  modified  LS  periodogram  method  for  each  of  these  12  sectors.  Knowing  the  phase 
distribution  of  the  wave  (the  variations  with  the  frequency  of  maximum  in  the  power  spestrum 
of  periodogram)  around  the  annulus,  the  propagation  direction  (azimuth)  and  horizontal  phase 
speed  can  be  determined  by  algorithm  described  by  Zhang  et  al.  (1993a).  This  algorithm  is 
based  on  the  last  squares  fitting.  Here  we  will  discuss  the  new  approach  to  the  analysis  of  the 
amplitude  and  phase  spectra  in  6  sectors  which  is  based  on  the  COPHASE  method  (Posmentier 
and  Herrmann,  1971;  Montes  and  Posmentier,  1971). 


COPHASE,  an  ad  hoc  statistic  characterizing  signal  strength  as  a  function  of  velocity  and 
direction,  for  data  from  an  array  of  detector  receiving  a  signal  and  uncorrelated  noise,  has  been 
defined  by  Posmentier  and  Herrmann  (1971)  as  follows 

E  E  E'[A,(a;„)  +  A,(a;„)]  cos[(^f  - 

n=7Vi  j=k-hl  k=l 

iV2  K  A'-l 

E  .E  E 

n=Ni  j=A;4*l  ^=1 

where  V  is  the  assumed  phase  speed;  0  is  the  assumed  azimuth  from  which  signal  arrives  (note 
that  in  our  realization  the  azimuth  indicate  the  direction  to  what  signal  propagates);  n  is  the 
index  of  frequency  (frequencies  in  the  band  from  the  Ai’th  to  the  A2’th  harmonic  are  used);  j 
is  the  station  number  (there  are  K  stations;  u;„  is  the  n’th  frequency;  Aj(a;„)  is  the  n’th  Fourier 
amplitude  at  the  y’th  detector;  is  the  n’th  Fourier  phase  at  the  y’th  detector;  <pj  is  the 
phase  of  the  n’th  frequency,  arriving  at  the  y’th  detector,  based  on  the  position  of  the  detector 
and  the  assumed  values  of  V  and  6. 


+  Ak{oJn)] 


Basically,  C{V,d)  calculates  the  normalized,  weighted  sum  over  all  \K{K  —  1)  detector  pairs, 
and  over  all  frequencies,  of  the  cosines  of  the  difference  between  the  Fourier  phase  differences 
and  the  assumed  phase  differences.  Note  that  C{V^9)  is  a  statistic  with  an  expectation  of  0 
if  the  K  records  are  non-correlated,  and  an  expectation  Of  -1-1.0  if  the  K  records  represent 
a  non-dispersive  propagating  signal,  at  the  assumed  value  of  V  and  6  corresponding  the  the 
actual  velocity  and  direction  of  the  signal.  COPHASE  technique  was  successfully  aplied  to  the 
analysis  of  seismic  Rayleigh  waves  and  acou'stic  signals  generated  by  the  Saturn-Apollo  rocket 
launches  (Montes  and  Postmentier,  1971).  In  the  case  of  MORTI  data  we  have  12  records  for 
each  sectors  and  COPHASE  technique  can  be  directly  applied  to  the  analysis  of  atmospheric 
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waves.  However,  to  exclude  the  errors  due  to  small  phase  differences  in  neighbour  sectors,  we 
used  only  6  pairs  of  opposite  sectors  to  calculate  C{V,9).  Besides,  in  the  case  if  the  amplitude 
of  considered  harmonic  fore  some  sector  is  small  in  comparison  with  the  amplitude  of  this  har¬ 
monic  in  the  record  of  data  averaged  over  12  sectors,  the  corresponding  pair  is  excluded  from 
the  consideration. 

The  relative  airglow  column  emission  rate,  O2  rotational  temperature,  and  results  of  LS  peri- 
odogram  analysis  applied  to  the  averaged  over  all  sectors  MORTI  data  obtained  at  Bear  Lake 
Observatory  on  January  20,  1993  are  shown  in  Fig.  9.  One  can  see  well  defined  maxima  in 
the  relative  emission  rate  and  temperature  power  spectra  at  the  period  of  180  min.  The  pe- 
riodogram  analysis  of  data  for  each  of  12  sectors  shows  the  presence  of  these  maxima  in  all 
time  series.  A  contour  plots  of  COPHASE  for  the  emission  rate  and  temperature  are  shown  in 
Fig.  10.  The  calculations  were  made  for  velocities  from  100  to  400  m/s  and  include  5  harmon¬ 
ics  (from  171  to  195  min  periods).  COPHASE  analysis  indicates  that  atmospheric  wave  was 
propagating  with  a  horizontal  phase  speed  of  about  250  m/s  and  that  azumuth  of  the  prop¬ 
agation  was  50°  ±  10°  counterclockwise  from  the  first  sector.  Unfortunately,  we  do  not  know 
now  what  was  the  orientation  of  the  MORTI  sectors  at  Bear  Lake  Observatory.  Fig.  10  also 
shows  that  the  results  obtained  by  COPHASE  analysis  applied  to  the  column  emission  rate  and 
temperature  data  are  very  similar.  The  values  of  azimuth  and  horizontal  phase  speed  obtained 
by  COPHASE  technique  and  last  squares  method  (Zhang  et  al.,  1993a)  are  presented  in  Table  3. 


Table  3 


Bear  Lake 

January  20-21,  1993 

Azimuth 

deg 

Phase  speed 
m/s 

Emission  rate 

COPHASE 

35 

250 

Last  squares 

33 

247 

Temperature 

COPHASE 

60 

250 

Last  squares 

25 

486 

This  table  shows  that  COPHASE  analysis  is  more  reliable  for  determination  of  the  azimuth  of 
propagation  and  horizontal  phase  speed  of  atmospheric  waves  then  last  squares  method  (the 
differences  between  the  values  obtained  from  the  emision  rate  and  temperatire  data  are  smaller). 
Using  the  value  of  the  horizontal  phase  speed  from  Table  3  we  can  estimate  the  horizontal  wave 
length  of  considered  disturbance  Xk  =  CftT=2700  km.  This  wave  length  is  significantly  larger 
then  the  distance  between  the  points  in  the  mesopause  region  from  what  we  registrate  the 
column  emission  rate  (for  opposite  sectors  this  distance  is  about  110  km).  In  result  the  phase 
shifts  are  very  small  (<  360°)  and  small  errors  in  determination  of  the  phase  due  to  influence 
of  long  period  harmonics  can  change  the  final  result  substantially.  To  avoid  the  influence  of 
long  period  harmonics  {T=  497.9,  288.3  min,  see  Fig.  9c),  we  removed  these  harmonics  from 
the  time  series  using  last  squares  fitting  and'xepeat  the  COPHASE  analysis.  The  relative  emis¬ 
sion  rate,  temperature,  and  results  of  LS  periodogram  analysis  applied  to  the  averaged  over  all 
sectors  data  after  removing  quadratic  trends  and  two  long  period  harmonics  from  time  series 
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Figure  9:  The  relative  emission  rate  (a)  and  rotational  temperature  (b).  The  normalized  LS 
periodograms  (c)  from  B'fBo  -  solid  line  and  T'^ITq  -  dashed  line  for  the  data  observed  at 
Bear  Lake  Observatory  on  Jan.  20  to  Jan.  21,  1993.  The  quadratic  trends  are  removed. 
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are  shown  in  Fig.  11.  The  results  of  COPHASE  analysis  applied  to  the  amplitude  and  phase 
spactra  for  12  sectors  time  series  after  similar  procedure  are  shown  in  Fig.  12.  The  values  of 
azimuth  and  horizontal  phase  speed  obtained  by  COPHASE  technique  and  last  squares  method 
(Zhang  et  ah,  1993a)  applied  to  the  data  after  removing  two  long  period  harmonics  from  time 
series  are  presented  in  Table  4. 


Table  4 


Bear  Lake 

January  20-21,  1993 

Azimuth 

deg 

Phase  speed 
m/s 

Emission  rate 

COPHASE 

45 

Last  squares 

43 

130 

Temperature 

COPHASE 

60 

165 

Last  squares 

44 

274 

The  comparison  of  Figures  10,  12  and  the  results  listed  in  Tables  3,  4  shows  that  long  period 
variations  strongly  influence  on  the  accuracy  of  the  phase  speed  determination  by  COPHASE 
and  last  squares  methods.  The  results  obtained  from  the  time  series  without  long  period  vari¬ 
ations  (Fig.  12,  Table  4)  show  that  true  values  of  horizontal  phase  speed  and  wave  length  are 
appoximately  half  as  much  as  the  values  obtained  from  initial  data.  Besides,  the  maxima  in 
COPHASE  analysis  become  better  defined  and  localized.  The  results  obtained  by  last  squares 
method  applied  to  the  emission  rate  time  series  show  good  agreement  with  COPHSE  analysis, 
but  the  determination  of  azimuth  and  phase  speed  is  incorrect  for  temperature  data.  Thus  we 
can  conclude  that  COPHASE  technique  is  more  reliable  and  useful  for  analysis  of  MORTI  data 
then  last  squares  method,  and  especially  this  may  be  important  if  we  have  not  so-good  time 
series  as  the  ones  presented  in  Figures  9,  11  (for  instance,  when  there  are  thin  tropospheric 
clouds  and/or  faint  background  illumination). 

Acknowledgements  -  We  are  grateful  to  Prof.  G.G. Shepherd  and  Dr.  R.H. Wiens  for  pro¬ 
viding  us  with  the  MORTI  instrument,  usefull  discussion  of  results,  and  for  sending  MORTI 
data  obtained  at  Bear  Lake  Observatory. 
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Figure  11:  As  in  Fig.  9,  but  two  long  period  harmonics  (T  =  497.9,288.3  min)  were  removed 
by  last  squares  fitting. 
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II.  lONOJSPHERIC  OBSERVATIONS 
2.1  The  BASIS  ionosonde  and  Doppler  sounder. 


Digital  multi-functional  ionosonde  called  BASIS  is  a  basic  instrument  which 
is  employed  for  ground  -  based  research  and  prediction  purposes  at  Institute  of 
Ionosphere.  BASIS  was  also  used  in  the  campaign  -  type  research  program  consid¬ 
ered  at  this  report.  BASIS  is  operated  as  a  vertical  incidence  sounder  giving  infor¬ 
mation  on  the  ionospheric  parameters.  BASIS  provided  direct  information  of  the  io¬ 
nospheric  signal  delays  (virtual  height  of  radio  pulse  reflections  called  ionograms)  as 
a  function  of  frequency.  The  main  parameters  of  the  ionosonde  are  presented  in  Ta¬ 
ble  5 


Table  5 


SPECIFICATIONS 

Pulse  Power 

15  kW 

Frequency  Range 

0.3  to  20  MHz 

Pulse  Coding 

Single  pulse 

Frequency  Sweeps 

400  frequencies 

Height  Range 

50  to  540  km 

Height  Resolution 

1.5  km 

BASIS  was  designed  to  record  ionograms  on  film  and  electrographic  paper. 
Limitations  connected  with  storing,  scaling  and  processing  of  ionogram  records 
made  us  to  begin  designing  a  computer  -  controlled  system  to  collect,  store,  display 
and  scale  the  ionograms  in  digital  form.  A  realization  of  these  fiinctions  was  made 
on  the  basis  of  a  "  Pentium  -  166"  computer.  The  ionosonde  and  the  computer  were 
coupled  by  a  printed  circuit  board  "  L  -  154".  L  -  154  passed  signals  from  the  re¬ 
ceiver  and  the  system  of  the  control  to  computer  in  order  to  collect  in  for  motion 
and  synchronize  operations  of  the  ionosonde  and  the  computer.  The  designed  sys¬ 
tem  allows  to  record  an  amplitude  of  reflected  signals  as  well.  In  order  to  improve  a 
quality  of  ionograms  the  changing  threshold  of  a  signal  detection  was  used.  Vari¬ 
ous  levels  of  the  threshold  allows  to  operate  under  different  levels  of  interference 
from  broadcasting.  A  software  package,  was  developed  to  implement  next  func¬ 
tions: 

a)  recording  of  signal  amplitudes  in  memory  of  the  computer 

b)  recording  of  ionograms  on  the  computer  screen  and  in  memory  of  the  computer. 
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c)  scaling  of  the  stored  data  by  arrow  keys. 

"  BASIS  "  sounds  4  times  at  each  frequency  .  Each  of  400  operating  frequen¬ 
cies  has  4096  levels  of  an  amplitude  quantization  and  360  levels  of  effective  height 
quantization.  The  signal  amplitudes  at  the  same  frequency  are  summed.  One  iono- 
gram  with  amplitudes  of  recorded  signals  requires  288  kB.  Each  ionogram  is  stored 
with  header,  date  and  time  information. 

Scaling  of  data  is  carried  out  by  using  a  semi  -  automatic  method.  O  and  X  - 
components  are  scaled  separately. 

Virtual  height  (h'(t))  temporal  variations  at  fixed  frequencies  can  be  obtained 
from  digital  ionogram.  These  variations  are  used  for  the  estimations  of  parameters  of 
traveling  ionospheric  distinbances.  To  obtain  the  h'(t)  temporal  variations  digital 
ionograms  were  recorded  every  5  minutes.  Fig.  13  gives  an  example  of  h'(t)  tempo¬ 
ral  Variations  at  several  fixed  frequencies.  These  variations  appear  to  be  caused  by 
the  passage  of  internal  atmospheric  gravity  waves  (AGWs)  through  the  ionosphere 
above  the  ionosoride.  The  variations  are  seen  to  travel  downward  indicating  well 
known  property  of  the  AGWs  vertical  phase  speed. 

Doppler  sounder  was  designed  in  Institute  of  Ionosphere  to  study  traveling 
ionospheric  disturbances  (Yakovets  et  al.,  1995).  The  equipment  uses  two  cross  - 
polarized  20  m  high  delta  type  antenna  for  transmission  and  reception.  A  frequency 
synthesizer  generates  four  predetermined  frequencies  (  f  i,  i=1...4)  which  were  se¬ 
lected  empirically  in  every  observation  in  order  to  ensure  an  approximately  equidis¬ 
tant  distribution  of  reflection  points  below  F  -  layer  peak  and  to  avoid  interference 
from  broadcasting  at  HF  frequencies. 

Reflected  from  the  ionosphere  signals  at  frequencies  f  i  +  f  di  (where  f  di  is  the 
Doppler  shift  caused  by  the  motion  of  the  ionospheric  screen)  were  received  and 
Doppler  shifts  were  computed  and  stored  for  fiirther  processing.  During  this  winter 
the  facility  was  equipped  by  two  channel  transformer  set  up  between  the  receiver 
and  computer  so  we  could  record  only  two  operating  frequencies.  Fig  .14  show 
Doppler  shift  variations  .  There  is  noticeable  difference  between  h'(t)  and  Doppler 
shift  variations.  The  first  ones  are  mostly  composed  of  long  period  variations  which 
are  typical  for  the  TIDs  night  observations  (Yakovets,  1995). 

Doppler  shift  variations  are  a  time  derivative  of  the  phase  path  of  radio  waves 
and  therefore,  Doppler  technique  enlarges  a  contribution  of  high  frequency  compo¬ 
nents  into  TIDs  spectra.  Simultaneous  observations  of  TIDs  by  the  ionosonde  and 
Doppler  soimder  showed  that  the  first  technique  was  more  informative  than  the 
second  one. 

So  the  ionosonde  was  considered  to  be  the  main  instrument  to  obtain  TIDs 
parameters  and  data  of  the  Doppler  sounder  played  an  additional  role. 
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2.2.  Experimental  results 

During  campaigns  conducted  from  October  1997  till  April  1998  the  iono¬ 
sphere  was  soimded  at  5  -  min  intervals.  Altogether,  over  5000  ionograms  were  re¬ 
corded  and  analysed.  Most  of  h'(t)  records  was  obtained  for  extra  ordinary  compo¬ 
nent  of  ionospheric  signals. 

In  order  to  study  a  process  of  TIDs  propagation  a  behaviour  of  virtual  heights 
(h'(t))  at  specific  frequencies  and  the  F  -  region  critical  frequencies  were  presented. 

He  lowest  specific  (operating)  frequency  was  chosen  to  be  of  2.5  MHz  be¬ 
cause  the  real  height  of  its  reflection  was  close  to  a  bottom  of  the  night  F2  -  region. 
The  adjacent  operating  frequencies  were  separated  by  0.5  MHz  in  October,  1997 
and  by  0.25  MHz  in  other  months. 

About  15%  of  the  total  number  of  h'(t)  records  were  rejected  from  the  analy¬ 
sis  because  of  a  presence  of  the  Eg  -  layer  which  screened  the  F  -  region  and  other 
reasons. 

Fig.  15  a  shows  typical  variations  of  h'(t)  on  the  set  of  operating  frequencies. 

In  our  Figures  we  use  local  time  which  leave  behind  75°  E  solar  mean  time  by 
1  hour  in  the  winter  time  (from  October,  27  till  March,  28)  and  2  hour  in  the  sum¬ 
mer  (from  March,  29  ).  Visual  examination  of  the  records  allows  to  do  some  con¬ 
clusions.  It  seems  obvious  the  existence  of  single  the  dominant  periodicity  with  a 
period  of  ~  2.2  hours  with  definite  TIDs  wave  structures.  Also  it  is  evident  of 
downward  phase  propagation  of  the  observed  waves  especially  at  high  operating 
frequencies  that  imply  upward  propagation  of  AGWs  energy  (Hines,  1960).  Addi¬ 
tionally  to  a  wave  process  observed  at  h'(t)  variations  there  is  a  slow  trend  which 
being  remained  in  data  can  distort  a  result  of  a  spectral  analysis. 

It  is  known  that  processes  in  ionospheric  plasma  and  neutral  atmosphere  drive 
the  height  of  the  F  -  region  (Bousanto,  1990,  Miller  et.  al.,  1986). 

After  sunset  median  heights  increase  and  continue  to  be  approximately  con¬ 
stant  during  the  night.  At  the  end  of  the  night  median  heights  decrease.  In  addition 
to  the  broad  night  -  time  plateau  in  the  h'(t)  running  the  morning  peak  occurs  fre¬ 
quently  close  to  the  F  -  region  sunrise  time.  Hajkowicz  (1991)  proposed  to  consider 
this  peak  as  a  result  of  AGWs  generation  by  the  supersonic  motion  of  the  sunrise 
terminator. 

In  sprite  of  an  apparent  regular  nature  of  the  night  h'(t)  trend  the  last  demon¬ 
strates  great  variety  of  its  behaviour.  Therefore,  we  could  not  approximate  it  by  a 
certain  fimction  in  order  to  remove  it  from  original  records.  A  filter  of  running 
means  was  used  to  remove  a  trend  from  the  row  records.  Choosing  a  length  of  the 
running  window  one  has  to  satisfy  to  two  contradictories.  From  the  one  hand  cut¬ 
ting  of  original  records  because  of  this  filtration  should  be  minumsed  and  from  the 
other  hand  the  length  of  the  running  window  is  recommended  to  be  enlarged  to  pro¬ 
vide  the  proper  filter  cut  -  off.  Visual  examination  all  of  row  h^t)  records  obtained 
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showed  that  most  of  period  of  dominant  waves  were  in  the  vicinity  of  2  hours  and 
this  value  was  chosen  to  be  a  length  of  the  running  window.  The  length  of  the  run¬ 
ning  window  is  considered  to  be  a  filter  cut-off  period,  but  data  processing  showed 
that  the  filter  cut-off  was  rather  smooth  and  therefore  even  when  periods  of  domi¬ 
nant  waves  exceeded  the  cut-off  period  they  preserved  a  significant  peak  in  spec¬ 
tral  estimations. 

Fig.  15  a  shows  that  an  amplitude  of  h'(t)  variations  grows  with  growing  of 
an  operating  frequency.  But  it  does  not  mean  that  an  amplitude  of  electron  density 
fluctuations  grows  in  the  same  way.  The  reason  of  this  is  following.  Commonly  the 
night  profile  of  an  electron  density  is  approximated  by  a  parabola  widi  is  a  good 
accuracy.  That  means  fast  rate  of  growing  of  the  electron  density  in  the  bottom  of 
the  F  -  region  and  slow  rate  near  the  maximum  of  the  F  -  region.  If  TIDs  propagate 
tlirough  the  ionosphere  with  relative  amplitudes  S  =  (No  -  electron  density,  AN 

-  ah  amplitude  of  electron  density  fluctuations  )  no  depending  of  a  height  this  cause 
small  variations  of  a  real  height  of  a  signal  reflection  in  the  bottom  of  the  layer  and 
large  ones  near  the  maximum  of  the  layer.  Those  fact  that  we  recorded  virtual 
height  variations  enhances  effect  of  apparent  growing  of  TIDs  amplitude  upward.  So 
virtual  height  variations  do  not  allow  to  define  a  TIDs  amplitude.  For  this  purpose 
critical  frequency  variations  were  used. 

Square  of  the  critical  frequency  of  the  ordinary  component  is  proportional  to 
the  electron  content  in  the  F  -  layer  maximum  (Budden,  1961): 


The  amplitude  of  electron  density  variations  ^  (half  of  the  peak-to-peak 
variations)  is  found  from  the  amplitude  of  critical  frequency  variations  as; 


N  /o 


(4) 


The  expression  (4)  can  be  used  to  estimate  the  relative  amplitude  of  electron 
content  variations  in  F  -  layer  maximum  caused  by  TIDs.  Fig.  15  b  shows  variations 
of  F  -  region  critical  frequencies  (foF,  fxF)  of  ordinary  and  extraordinary  components 
of  the  ionospheric  signal  for  the  same  night  as  in  Fig.  15  a.  Comparing  fr(t)  and  foF, 
fxF  variations  one  can  see  that  the  last  are  about  in  antiphase  with  h'(t)  ones. 

This  effect  is  always  present  in  the  experimental  records.  It  seems  can  be  ex¬ 
plained  by  a  large  difference  between  estimated  vertical  wavelength  (X.z)  and  a  half 
thickness  of  the  night  F-  layer.  An  average  distance  between  a  height  of  the  night  F 
-  layer  maximum  and  its  bottom  does  not  exceed  ~50  km.  Taking  into  account  that 
average  ?iz=250  km  (Table  6)  one  can  conclude  that  the  electron  density  below  the 
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F  -  layer  ma^mum  fluctuates  practically  synchronously  at  all  heights.  That  means 
that,  for  example,  decreasing  of  the  electron  density  in  the  F  -  layer  maximum 
caused  by  wave  propagation  accompanies  with  electron  decreasing  at  all  F  -  layer 
heights.  This  decreasing  shifts  points  of  a  reflection  for  all  operating  frequencies 
upward  and,  as  a  result,  antiphase  variations  of  fo,  x  F(t)  and  h'(t)  are  observed. 

In  order  to  obtain  power  spectra  and  phase  differences  between  h'(t)  varia¬ 
tions  occurred  at  various  height  the  method  of  Blackman  and  Tukey  (lenkins  and 
Watts,  1968)  has  been  employed.  The  Tukey  window  was  used  as  a  lag  window. 
The  maximum  lag  of  120  min  has  been  adopted  to  obtain  a  reasonable  frequency 
resolution  and  accuracy  of  phase  differency  estimations.  Fig.  16  shows  power  spec¬ 
tra  for  h'(t)  records  of  operating  frequencies  ft  and  fs  shown  in  Fig.  15a. 

The  ft  is  the  lowest  operating  frequency  and  the  fs  is  the  fifth  operating  fre¬ 
quency  spaced  by  16  km  from  the  ft.  Power  spectra  show  a  prominent  peaks  at  the 
frequency  Fs  0.54  cph  (cycles  per  our)  and  their  positions  on  the  frequency  axis 
are  very  close. 

These  spectra  are  typical  for  h'(t)  records  obtained.  They  consist  of  a  single 
prominent  peak  corresponding  to  dominant  wave.  Small  peaks  on  the  lower  periods 
are  often  present  but  their  values  are  significantly  less. 

Not  all  of  h'(t)  power  spectra  demonstrates  a  coincidence  of  spectral  peaks. 
Fig.  17  shows  spectra  h'(t)  variations  at  frequencies  ft  and  fs  obtained  in  December 
25,  1987.  It  is  seen  that  a  spectral  peak  at  the  lower  (ft)  operating  frequency  is 
slightly  shifted  toward  a  shorter  period. 

T.Shibata  and  K.  Schlegel  (1993)  also  observed  a  frequency  shift  of  spectral 
peaks  and  gave  a  reasonable  explanation  of  this  fact.  TIDs  propagate  on  the  back¬ 
ground  of  a  neutral  horizontal  wind  which  shift  an  TIDs  intrinsic  frequency  due  to  a 
Doppler  effect.  The  Doppler  effect  of  a  horizontal  wind  may  be  expressed  in  the 
following  fonn: 


f  =  F  +  —k’U{K) 

In 

Here  F'is  a  TID  apparent  frequency  which  is  recorded  by  an  observer  on  the 
ground,  F  is  a  TID  intrinsic  frequency,  k  is  wave  number  vector  and  u{h)  is  a 
horizontal  wind  velocity  vector  Assuming  that  the  intrinsic  wave  parameters,  F 
and  k  ,  are  independent  of  a  height,  one  can  conclude  that  a  vertical  gradient  of  the 
horizontal  wind  velocity  vector  causes  a  height  dependence  of  apparent  frequency. 
It  is  seen  from  (4)  that  a  set  of  F  for  series  heights  gives  a  possibility  to  compute  a 
vertical  gradient  of  the  horizontal  wind  if  k  is  known.  We  could  not  define  k  be¬ 
cause  we  sounded  the  ionosphere  from  the  single  point. 
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Frequency  shifts  between  positions  of  spectral  peaks  have  to  be  accompanied 
with  gradual  temporal  changing  the  phase  difference  between  h'(t)  variations  oc¬ 
curred  at  different  operating  frequencies.  Fig  .18  gives  a  set  of  h'(t)  variations  which 
spectra  are  plotted  in  Fig.  17  Visual  check  confirms  this  fact.  One  can  see  that  a 
phase  lag  between  the  same  h^t)  signatures  at  various  operating  frequencies 
changes  in  time.  Vary  fast  changing  takes  place  between  tlie  first  and  second  crests 
of  the  waves.  Next  crests  reveal  more  moderate  changing  in  the  phase  lag.  These 
noticeable  fluctuations  of  phase  lag  changing  are  smoothed  by  a  technique  of  proc¬ 
essing  which  give  an  average  rate  of  phase  lag  changing.  Temporal  fluctuations  of 
phase  lag  changing  may  be  caused  by  fluctuations  of  a  vertical  gradient  of  the  hori¬ 
zontal  wind. 

From  the  study  of  data  obtained  it  is  clear  that  TIDs  are  practically  always 
present  in  the  night  ionosphere.  Now  we  will  classify  TIDs  according  to  their  be¬ 
havior  during  the  period  of  observations.  The  first  type  presents  TIDs  existing  dur¬ 
ing  all  night  with  approximately  constant  amplitude  and  period  (Fig.  13,  15a).In 
these  figures  gaps  iii  the  run  of  the  top  operating  frequencies  correspond  to  periods 
when  the  critical  frequency  becomes  below  the  operating  frequency. 

The  second  type  presents  TIDs  with  large  amplitudes  in  the  beginning  of  rec¬ 
ords  and  their  further  gradual  attenuation  (Fig.  19,20).  These  two  types  of  h’(t) 
variations  are  the  most  frequently  met  in  the  analyzed  observations.  This  types  ap¬ 
pear  to  be  understood  on  the  basis  of  packet-like  TIDs  structure  considered  in  sec¬ 
tion  2.3. 

Kaliev  et  al.  (1989)  analyzing  packet  with  the  periods  of  oscillations  of  7-60 
min.  found  a  correlation  between  a  period  of  oscillations  and  a  duration  of  a  packet. 
Long  period  packets  could  last  as  long  as  6  hours.  Hence,  a  duration  of  packets  with 
a  period  of  oscillations  of  two  hours  may  reach  a  significant  value  exceeding  a 
course  of  observations.  We  can  think  that  nights  with  the  constant  amplitude  of  h’(t) 
variations  correspond  to  a  middle  fraction  of  a  wave  packet.  The  second  type  of  the 
amplitude  variations  behavior  may  present  the  end  of  current  packet. 

Fig.21  shows  small  variations  of  h’(t)  during  the  night  and  large  surge  in  the 
premomi^  ^urs.  This  surge  seems  to  be  caused  by  the  solar  terminator  generating 

near  sunrise  on  the  height  of  the  F-region. 

In  order  to  summarize  TIDs  characteristics  measured  during  the  period  of  ob¬ 
servations  a  table  6  was  made. 
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25.12.97 

96 

30.12.97 

90 

31.12.97 

180 

5.01.98 

108 

23.01.98 

102 

29.01.98 

126 

30.01.98 

108 

368 

no 

122 

+ 

171 

no 

2.02.98 


24.02.98 


27.02.98 


28.02.98 


2.03.98 


23.03.98 

24.03.98 


29.03.98 


30.03.98 


2.04.98 


no 

no 

+ 

■ 

no 

no 
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The  first  column  contents  dates  of  observations  analyzed.  The  second  one 
gives  a  period  of  a  dominant  wave  observed  in  the  certain  night.  The  sign  ‘’uncert.” 
Denotes  uncertain  data.  Relative  amplitudes  computed  from  formula  (4)  are  listed  at 
the  third  column.  Here,  the  sign  ‘’uncert.”  Corresponds  to  events  when  wavC  os¬ 
cillations  occurred  on  the  level  or  below  of  a  noise  amplitude  which  average  value 
was  estimated  to  be  of  2  %.  The  fourth  column  contents  vertical  phase  speeds  esti¬ 
mated  from  a  difference  of  real  heights  of  ionospheric  signal  reflections  and  lags 
obtained  from  cross-correlation  analysis.  For  seven  nights  speeds  could  not  be  esti¬ 
mated  because  very  low  time  delays.  Four  of  these  events  corresponded  to  low 
wave  amplitude  may  be  caused  by  errors  of  statistical  processing.  But  three  events 
(29.03.98,  30.03.98,  2.04.98)  for  wave  amplitudes  exceeding  noise  levels.  Here, 
low  lags  could  be  explained  on  the  basis  of  Perkins  (1973)  instability.  Perkins  insta¬ 
bility  develops  in  the  night  F-region  and  force  the  ionosphere  to  move  up  and  down 
as  whole..  So  the  vertical  phase  speed  estimated  from  record  obtained  for  different 
operating  frequencies  would  be  uncertainly  large..  Kelley  and  Fukao  (1991)  pro¬ 
posed  that  gravity  waves  in  the  certain  conditions  can  initiate  the  Perkins  instability. 

Next  column  gives  values  of  the  vertical  wavelength  estimated  from  the  prod¬ 
uct  of  a  period  by  Vz . 

The  column  with  a  frequency  shift  gives  a  difference  between  the  spectral 
peaks  at  operating  frequencies  fi  and  fs .  The  positive  sign  corresponds  to  an  event 
when  a  spectral  peak  at  the  lower  (fi)  operating  frequency  is  shifted  toward  a 
shorter  period.  The  negative  sign  corresponds  to  the  opposite  direction  of  a  fre¬ 
quency  shift.  And  sign  ‘’no”  corresponds  to  events  when  no  shifts  were  observed. 

In  the  last  column  Kp-indices  of  geomagnetic  activity  are  listed. 

Parameters  from  the  Table  6  are  presented  by  distributions  in  the  Fig.22. 
Plots  show  that  the  main  part  of  parameters  analyzed  is  distributed  in  relatively  nar¬ 
row  ranges.  Fig.  23  presents  a  relative  amplitude  of  TIDs  (8)  versus  Kp-index.  There 
is  not  evident  correlation  between  these  values,  but  one  can  see  that  under  high 
geomagnetic  activity  wave  amplitudes  noticeably  exceeded  the  noise  level.  This 
confirm  results  of  TIDs  observations  after  great  geomagnetic  storm  (Hajko- 
wicz,I987,  1991,  Yeh  et  al.,  1994).  At  the  same  time  for  nights  with  low  geomag¬ 
netic  activity  wave  amplitudes  are  distributed  in  the  rather  wide  range.  Large  am¬ 
plitudes  for  the  low  Kp-indices  prove  that  there  are  sources  of  TIDs  activity  not 
connected  with  auroral  regions. 

Fig.  23b  presents  a  seasonal  dependence  of  an  average  relative  amplitude.  It 
grows  noticeable  in  the  winter  months. 
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To  compare  the  results  of  measurements  in  the  ionosphere  and  mesosphere 
three  days  were  chosen  with  wave  motions  well  developed  into  both  regions.  De¬ 
trending  applied  for  processing  of  optical  data  was  used  for  analysis  of  the  iono¬ 
spheric  data.  This  procedure  allowed  to  reveal  the  low  frequency  components  elimi¬ 
nated  by  cut-off  high  frequency  filtering.  Fig.  24,  25,  26  give  detrended  relative  vir¬ 
tual  height  variations  at  frequencies  fl,  f5  and  their  normalized  power  spectra  for 
Nov.  27,  1997,  Feb.  2,  1998  and  March  23,  1998.  Here,  a  solid  line  corresponds  to 
fl  and  a  dashed  line  -  to  f5.  Temporal  variations  of  relative  emission  rate  tempera¬ 
ture  and  their  spectra  are  plotted  in  Fig.  2,  6^  8.  Comparing  ionospheric  and  meso¬ 
spheric  wave  variations  one  can  found  their  similarity  and  this  fact  is  supported  by 
their  spectra  except  the  spectra  for  Nov.  27,  1997.  The  optical  spectra  show  the 
very  large  peaks  at  period  of  500  min.,  but  the  ionospheric  spectra  have  no  the 
similar  peaks.  This  fact  is  explained  by  different  lengths  of  optical  and  ionospheric 
data.  The  length  of  optical  data  equaled  to  12  hours  allows  to  resolve  the  500  min 
component.  Meanwhile  the  length  of  ionospheric  data  for  this  date  is  equal  to  8 
hours,  therefore  the  500  min  component  was  considered  as  trend  and  was  eliminated 
from  row  data.  In  order  to  equalize  the  pre-processing  procedure  the  same  detrend¬ 
ing  was  applied  to  optical  data  and  corresponding  data  are  plotted  in  Fig.27.  Peri¬ 
ods  of  prominent  spectral  peaks  for  virtual  height  (h’f5,  h’fl)  and  emission  rate  (E) 
variations  plotted  in  Fig.  24-27  are  listed  in  Table  7. 

Table  7 


Nov.27,  97 

Feb.  2,  98 

March  23,  98 

Tl(min) 

T2(mm) 

Tl(niin) 

T2(min) 

Tl(min) 

T2(min) 

h’f5 

143 

300 

110 

195 

130 

310 

h’fl 

135 

265 

no 

195 

120 

360 

E 

160 

300 

100 

150 

140 

300 

One  can  see  that  there  is  no  the  complete  coincidence  for  periods  of  spectral  peaks 
corresponding  to  different  heights  of  observations  except  ionospheric  data  for  Feb. 
2,  1988.  We  suppose  that  this  difference  is  caused  by  a  vertical  gradient  of  the 
background  neutral  wind,  which  shifts  the  intrinsic  frequency  for  different  heights 
on  a  different  value. 

The  comparison  of  optical  and  ionospheric  data  proves  simultaneous  existence  of 
wave  variations  with  close  periods  at  mesospheric  and  ionospheric  heights.  The  di¬ 
rection  of  the  vertical  phase  speed  of  a  wave  motion  with  a  middle  period  of  2  hours 
at  ionospheric  heights  (the  average  real  height  of  ionospheric  observations  was 
equal  to  250  km)  showed  that  these  are  travelling  ionospheric  disturbances  caused 
by  AGW’s.  It  is  of  great  interest  for  determine  a  nature  of  ionospheric  and  meso- 
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spheric  variations  with  periods  of  ~  300  min  whether  they  belong  to  AGW’s  of 
tidal  waves. 

Conclusions. 

Observations  of  travelling  ionospheric  disturbances  during  tlie  period  from 
October  1997  till  March  1998  showed  that  TIDs  with  periods  of  80-180  min  are 
common  feature  of  the  night  ionosphere.  Various  parameters  of  TIDs  including  pe¬ 
riods,  relative  amplitudes,  vertical  phase  speeds,  vertical  wavelengths  were  meas¬ 
ured  and  their  distributions  were  constructed.  A  frequency  shift  between  TID  spec¬ 
tral  peaks  belonging  to  different  heights  was  found  for  a  part  of  observations  and  its 
possible  reason  was  discussed.  A  seasonal  dependence  of  TIDs  average  relative 
amplitudes  for  six  months  was  plotted.  Evident  correlation  between  TID  relative 
amplitudes  and  Kp-index  was  not  found.  A  comparison  ionospheric  and  meso¬ 
spheric  data  showed  that  waves  with  nearby  frequencies  exist  simultaneously  in  the 
both  regions. 

Several  reasons  make  us  to  proceed  ionospheric  and  mesospheric  observations 
till  October  1988  to  complete  a  year  cycle.  We  hope  to  obtain  a  seasonal  depend¬ 
ence  of  various  waves  parameters.  Besides,  summer  time  is  favorable  for  studying 
travelling  planetary  waves  in  the  mesopause  region,  and  we  are  going  to  perform 
regular  optical  measurements  this  summer  if  even  the  period  of  observations  during 
the  night  will  be  short. 

Our  measurements  will  be  supported  by  the  optical  measurements  at  other  sites 
(Canada,  Spain,  Japan)  and  we  will  have  a  good  possibility  for  analysis  of  planetary 
waves.  Moreover,  several  tasks  connected  with  computations  have  to  be  done.  It  is 
necessary  to  compute  the  vertical  phase  speed  of  low  frequency  wave  components 
(T  ~  300  min.)  and  define  their  nature.  Calculations  of  the  wave  horizontal  velocity 
on  the  basis  of  optical  data  have  to  be  done  and  vertical  gradient  of  the  neutral  wind 
will  be  obtained. 

If  the  night  TIDs  spectra  reveal  a  prominent  peaks  corresponded  to  dominant  waves 
then  day-time  spectra  seem  to  be  broad  and  extended  toward  high  frequencies.  This 
fact  gives  a  possibility  to  study  a  fine  structure  of  TIDs  and  relationships  between 
different  parts  of  the  spectra.  The  results  of  this  study  will  be  present  in  sections  2.3, 

2.3.  A  STUDY  OF  WAVE  PACKETS  IN  TRAVELLING  IONOSPHERIC  DISTERBACES 

Hines  (1960)  was  the  first  who  suggested  that  internal  atmospheric  gravity  waves 
(AGWs)  propagating  in  the  neutral  thermospheric  gas  caused  travelling  ionospheric 
disturbances  (TIDs).  Many  radio  techniques  including  ground-based  and  satellite 
ones  were  employed  to  observe  TIDs  during  day-  and  night-time  in  the  different 
seasons  and  under  various  geomagnetic  and  solar  conditions. 
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There  were  at  least  two  approaches  to  an  interpretation  of  TIDs  observations. 
According  the  first  one,  only  TIDs  visually  identified  with  monochromatic  individ¬ 
ual  AGWs  were  taken  into  consideration  in  order  to  estimate  their  parameters 
(Davies  and  Jones,  1971;  Waldock  and  Jones,  1987).  This  approach  gave  an  op¬ 
portunity  to  obtain  TIDs  characteristics  fi'om  a  small  fi^action  of  experimental  data. 
For  example,  Jacobson  and  Carlos  (1991)  could  find  only  a  few  monochro¬ 
matic  TIDs  fi'om  their  database  which  exceeded  100  hours  of  observations. 

According  the  second  approach,  TIDs  are  considered  to  be  a  stationaiy  field 
with  certain  spatial  and  temporal  spectra.  In  order  to  estimate  TIDs  parameters  at 
all  fi’equencies  of  spectra  the  cross-spectral  analysis  was  employed  (Herron,  1973; 
Tsutsui  et  al,  1984;  Shibata  and  Okuzava,  1983,  Yakovets  et  al.,  1995).  This 
technique  allowed  to  utilise  practically  all  Doppler  and  other  continuous  TIDs  rec¬ 
ords  but  was  not  able  to  separate  monochromatic  individual  TIDs  and  background 
ones  which  were  considered  to  be  a  result  of  a  multiple-wave  overlap  that  distort  the 
spatial  and  temporal  coherence  (Jacobson  and  Carlos,  1991). 

Passing  TIDs  records  through  pass-band  filters,  Beggs  and  Butcher  (1989) 
and  Kaliev  et  al.  (1989)  have  found  that  some  fi’equency  ranges  had  temporal  sec¬ 
tions  which  comprised  packet-like  structures  with  amplitudes  of  oscillations  that 
appreciably  exceeded  the  amplitudes  of  background  oscillations.  JACOBSON  et  al. 
(1995)  offered  to  use  temporal  sections  with  large  amplitudes  to  estimate  the  wave 
parameters  because  they  are  thought  to  be  caused  by  individual  TIDs,  while  the 
background  oscillations  were  a  result  of  a  multiple  TIDs  superposition.  They  looked 
the  records  of  radio  signals  from  geosynchronous  satellite  beacons  through  time- 
and  frequency-windows  in  order  to  select  dominated  single-wave  individual  TIDs 
and  define  azimuths  of  their  propagation. 

It  should  be  pointed  out  that  the  time-  and  frequency-window  procedure  al¬ 
lows  to  find  not  only  parameters  of  individual  TIDs,  but  parameters  of  multiple 
TIDs  superposition  as  well.  A  comparison  of  parameters  of  individual  and  multiple 
TIDs  is  of  great  interest  because  it  gives  a  tool  to  inspect  previous  measurements 
interpreted  in  the  framework  of  a  model  of  stationary  TIDs. 

This  paper  presents  results  of  this  comparison.  The  technique  of  a  computa¬ 
tion  of  the  cross-correlation  coefficient  for  individual  TIDs  existing  on  the  back¬ 
ground  of  a  multiple  TIDs  superposition  was  developed.  This  technique  was  em¬ 
ployed  to  TIDs  observations  obtained  by  the  multifrequency  HF  Doppler  sounder. 
Besides  the  vertical  phase  speed  for  individual  TIDs  was  estimated  in  comparison 
with  the  speed  of  the  TIDs  superposition  and  a  behaviour  of  a  current  phase  was 
studied. 

Observations  of  TIDs  at  several  heights  of  the  F-layer  ionosphere  were  carried 
out  at  the  Alma-Ata  multifrequency  HF  Doppler  sounder  (76®55'E,  43®15'N).  The 
sounder  operated  at  five  sounding  frequencies  which  were  selected  for  every  obser¬ 
vation  to  ensme  an  approximately  equidistant  vertical  distribution  of  reflection 
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points  below  the  F-layer  peak.  The  details  of  the  equipment  have  been  described  in 
our  previous  paper  (Yakovets  at  al,  1995).  Fig.2  of  this  paper  shows  typical  rec¬ 
ords  of  the  Doppler  variations  (Trf  (0  )  four  operating  frequencies.  In  order  to  ob¬ 
tain  the  power  spectra  the  method  of  Blackman  and  Tukey  (JENKINS  and  WATTS, 
1968)  has  been  employed.  The  Tukey  window  was  used  as  a  lag  window.  The 
maximum  lag  of  60  min  has  been  adopted  to  obtain  a  reasonable  frequency  resolu¬ 
tion.  Figure  3a  (Yakovets  at  al.,  1995)  shows  normalised  power  spectra  for  the  re¬ 
cords  of  Doppler  variations  presented  in  Fig.  2a.  Careful  inspection  of  the  normal¬ 
ised  power  spectra  shows  that  typical  spectra  often  comprised  several  spectral 
peaks  on  the  background  of  continuous  spectra. 

In  order  to  study  the  temporal  behaviour  of  spectral  contributions  whose 
amplitudes  exceeded  the  background  ones,  the  complex  demodulation  method 
(Marple,  1988)  was  employed  .  This  method  operates  as  a  band-pass  filter  giving 
an  opportunity  to  obtain  both  a  filtered  signal  and  a  current  phase  of  this  signal  on 
its  output.  The  current  phase  is  present  a  difference  between  a  phase  of  the  filtered 
signal  and  a  reference  signal  with  a  fixed  frequency  equalled  to  the  central  fre¬ 
quency  of  the  filter.  The  filter  band-width  of  1  cph  (cycle  per  hour)  was  chosen  to 
be  close  to  spectral  resolution  of  the  method  of  Blackman  and  Tukey  used  at  the 
preliminary  stage  of  processing  and  to  ensure  a  good  reconstruction  of  wave  packets 
whose  presence  into  original  records  was  not  usually  apparent  to  the  eye  especially 
in  the  high  frequency  TDDs  range.  Figure  28  is  the  result  of  filtering  of  Doppler  shift 
variations  at  five  operating  frequencies  for  the  day  of  17  March  1988.  One  can 
make  the  certain  conclusions  about  a  temporal  behaviour  of  filtered  records. 

In  the  first  time  interval  {t  =  09.00  12.20  LT)  of  plotting  fd(t)  variations 

there  is  a  narrow  band  quasi-stochastic  process  and  in  the  second  time  interval  (/  = 
12.20  14.20  LT)  fci(t)  oscillations  have  the  appearance  of  a  wave  packet.  By  using 
this  technique  of  processing  we  found  that  packet-like  stmcture  in  fjt)  variations 
appeared  to  be  common.  It  is  therefore  important  to  compare  features  of  these  two 
types  of  TIDs  separately.  A  term  , ’’quasi-stochastic  TIDs”,  is  used  instead  of  ”a 
multiple  TIDs  superposition”  to  reflect  a  random  character  of  a  behaviour  of  an 
amplitude  and  phase  of  the  wave  resulted  from  the  multiple  TIDs  superposition. 

The  cross-correlation  fimctions  (  R{t)  )  between  filtered  records  of  fd{t)  were 
computed  for  quasi-stochastic  parts  of  the  records  and  for  packet-like  structures 
separately.  Fig.  29  presents  cross-correlation  functions  the  first, 

quasi-  stochastic  parts  of  records  plotted  in  Fig.  28.  Fig.  30  presents  cross¬ 
correlation  functions  the  second,  packed-like  parts  of  the  same  rec¬ 
ords.  From  these  figures  it  can  be  seen  that  values  of  peaks  of  exceed 

those  of  the  correspondent  peaks  of  Besides  lags  of  peaks  rela- 
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tive  to  r  =  0  exceed  lags  of  correspondent  „^(r)  peaks.  That  means  that  a  verti¬ 
cal  phase  velocity  of  packets  is  less  than  a  velocity  of  quasi-stochastic  TIDs. 

Values  of  the  vertical  phase  velocity  were  computed  from  the  measured 
heights  of  radiowaves  reflection  points  and  lags  of  correspondent  R{t)  peaks.  The 
true  heights  of  the  reflection  points  for  each  observation  were  computed  by  using 
ionograms  of  the  Alma-Ata  ionospheric  station.  All  registrated  waves  had  a 
downward  phase  motion,  and  their  velocities  grew  with  growing  of  a  height  of  the 
radio  wave  reflection.  In  Table  8  average  velocities  computed  from  records  corre¬ 
sponding  to  maximum  and  minimum  heights  for  dates  indicated  are  listed. 

One  can  see  that  for  all  considered  records  the  vertical  phase  velocity  of 
quasi-stochastic  TIDs  markedly  exceeded  of  tlie  velocity  of  packets.  Table  8  con¬ 
tents  velocities  for  both  packets  and  quasi-stochastic  parts  only  for  long-period 
TIDs.  The  reason  of  this  will  be  discussed  later. 

The  value  of  the  cross-correlation  function  peak  is  the  cross-correlation  coef¬ 
ficient.  It  gives  a  measure  of  the  coherence  between  records  corresponding  to  dif¬ 
ferent  ionospheric  heights.  It  is  clear  that  the  cross-correlation  coefficient  computed 
for  time  intervals  characterising  by  a  packet-like  structure  depends  not  only  on  the 
degree  of  the  coherence  of  wave  packets  but  and  a  ratio  of  amplitudes  of  packets 
and  quasi-stochastic  TIDs. 

Indeed,  caused  by  numerous  soiuces  the  quasi-stochastic  TIDs  field  insures 
continuous  variations  of  measured  7^/)  at  any  time  and  place  of  observations,  thus 
wave  packets  appear  only  on  the  backgroimd  of  quasi-stochastic  variations.  In  order 
to  obtain  the  cross-correlation  coefficient  i^^^  ^^for  the  “pure”  wave  packets,  the 

method  of  the  definition  of  Rs^  s^  was  developed  (Appendix).  By  using  this  method 
Rs^  Sj  were  estimated  and  listed  in  Table  8.  One  can  see  that  the  less  correlated 
quasi-stochastic  TIDs  behave  as  an  unwanted  noise  which  cause  a  drop  of  the  re¬ 
sultant  cross-correlation  coefficient  .  The  results  presented  may  be  di¬ 

vided  into  two  groups.  For  the  first  group  comprising  long  period  TIDs  the  values 
of  R^^  ^^  and  R„^  ^  are  sufficiently  large,  and  Rs^  s^  exceeds  them  insignificantly. 

The  different  situation  one  finds  for  the  second  group  comprising  short  period  TIDs 
(  T  <  20  min).  Drastic  increasing  of  R^^  s^  relative  to  and  especially  R„^  „^ 

occurred.  That  means  that  quasi-  stochastic  TIDs  mask  the  fact  that  even  high  fre¬ 
quency  TIDs  packets  propagate  practically  without  changing  their  form  along  the 
large  distance  on  the  vertical  plane. 

We  could  not  obtain  accurate  values  of  the  vertical  velocities  for  the  TIDs  pe¬ 
riod  less  then  20  minutes  from  cross-correlation  functions  because  of  small  lags  of 
their  maximums  relative  to  0  which  were  of  the  order  of  the  sample  interval  of  1 
min.  In  order  to  estimate  vertical  phase  speeds,  current  phases  were  computed  by 
the  complex  demodulation  method.  An  example  of  filtering  data  obtained  on  14  Feb. 
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1990  is  plotted  in  Fig.  31a.  The  current  phase  of  figured  data  are  plotted  in  27b.  The 
current  phases  corresponding  to  adjacent  heights  are  spaced  2  n  apart  to  avoid  a 
mishmash  of  their  runs. 

Comparing  Fig.  31a,  31b  one  can  see  two  successive  wave  packets  on  tempo¬ 
ral  intervals  12.20  13.40  and  13.50  -r  14.50.  There  is  a  significant  difference  be¬ 

tween  a  phase  behaviour  inside  and  outside  wave  packets.  If  a  current  phase  inside 
wave  packets  keeps  a  stable  position  with  small  fluctuation,  then  outside  the  phases 
fluctuate  with  a  large  amplitudes  and  independently  for  different  sounding  fi-equen- 
cies.  Small  fluctuations  inside  the  wave  packets  appear  to  be  caused  by  the  addi¬ 
tional  quasi-stochastic  oscillations.  Sometimes  a  gradual  run  of  the  current  phase  is 
distorted  by  its  fast  jumps.  Two  types  of  the  phase  jumps  were  found.  The  first  one 
is  characterised  by  jumps  which  are  less  of  In,  and  their  nature  will  be  discussed 
further.  The  second  one  is  characterised  by  jumps  equalled  to  2%,  and  they  originate 
from  the  phase  uncertainty  of  27c.  A  number  of  these  jumps  can  be  decreased  by 
certain  changing  of  the  initial  phase  of  the  records  to  shift  the  most  probable  phase 
in  the  filtered  record  in  the  vicinity  of  n,  but  it  is  difficult  to  eliminate  the  jumps 
completely  because  of  large  phase  fluctuations  into  quasi-stochastic  parts  of  the 
records. 


After  averaging  the  current  phase  inside  packets,  one  can  find  an  average  phase 
difference  between  various  heights,  and  estimate  the  vertical  phase  velocity.  This 
procedure  can  not  be  applied  to  temporal  intervals  with  quasi-stochastic  TIDs  be¬ 
cause  of  large  uncorrelated  phase  fluctuations  during  interval  09.00  ^  12.20  and 
14.50  16.30.  Therefore  for  TIDs  within  period  less  than  20  min  the  vertical  phase 

velocity  was  estimated  only  for  the  wave  packets. 

Fig.  32a,b  show  filtered  records  and  current  phases  at  the  frequency  of  3  cph 
on  14  February  1990.  Data  presented  in  Fig.  32  are  obtained  from  the  same  records 
as  in  Fig.  31.  However,  their  character  is  markedly  different.  Practically  the  whole 
length  of  records  is  occupied  with  three  distinct  wave  packets.  Looking  at  the  slope 
of  the  current  phase,  one  can  make  several  conclusions: 

a)  the  frequencies  of  oscillations  in  all  packets  are  slightly  different  from  the  central 
frequency  of  the  filter.  It  can  be  seen  from  the  behaviour  of  the  derivative  of  the 
phase  with  respect  to  time  that  defines  a  difference  between  the  central  fre¬ 
quency  of  the  filter  and  a  frequency  of  the  filtered  signal; 

b)  the  frequency  of  oscillations  of  the  second  and  third  packets  are  equal;  moreover 
they  are  coherent  because  their  frequencies  are  stable  in  the  course  of  a  packets 
existence; 

c)  clear  discontinuities  occur  between  the  second  and  the  third  packets,  and  a  value 
of  the  phase  jump  grows  witli  the  height  of  reflection  of  radio-signals; 


|j[^9  frequency  of  oscillations  of  the  first  packet  does  not  coincide  with  the  fre- 
^bncy  of  oscillations  of  the  second  packet; 
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e)  phase  differences  between  packets  corresponded  to  different  heights  are  reasona¬ 
bly  constant  except  frequency  fi ,  where  relatively  small  packet  amplitude  is  dis¬ 
torted  by  an  additional  noise. 

Table  8 


N 

Date 

71[min) 

F(ms-‘) 

Cross-  correl.  coefficient 

pac¬ 

ket 

stoch. 

field 

Xi,X2 

^Sl  yS2 

Ml 

km. 

1. 

17.03.88 

40 

100 

195 

0.94 

0.98 

0.83 

35 

2. 

08.02.90 

30 

175 

220 

0.95 

0.99 

0.85 

55 

3. 

14.02.90 

50 

165 

390 

0.87 

0.96 

0.70 

70 

wm 

19.11.90 

.40 

105 

165 

0.99 

0.99 

0.82 

50 

5. 

24.02.90 

18 

470 

0.80 

0.99 

0.4 

45 

6. 

19.11.90 

13 

350 

0.85 

0.98 

0.6 

50 

7. 

17.03.88 

7 

580 

0.78 

0.83 

0.7 

35 

8. 

14.02.90 

12 

500 

0.8 

0.98 

0.4 

70 

9. 

08.02.90 

11 

660 

0.6 

0.78 

0.35 

55 

10. 

15.02.90 

20 

650 

0.78 

0.99 

0.4 

60 

The  properties  of  the  individual  wave  packet  phase  structure  are  important 
because  they  may  help  to  realise  both  the  wave  generation  mechanisms  and  condi¬ 
tions  of  the  wave  propagation.  We  studied  28  series  of  filtered  daytime  records 
comprising  49  series  of  wave  packets.  Every  series  of  the  filtered  records  con¬ 
sisted  of  4-5  records  corresponding  to  different  operating  radio  frequencies  .  The 
average  length  of  the  records  was  5-6  hours.  Under  the  series  of  wave  packets  we 
meant  a  wave  packet  observed  at  all  heights  simultaneously. 

Most  wave  packets  had  a  constant  period  of  oscillations,  and  only  five 
packet  series  showed  a  growing  period  with  time  as  Francis  (1974)  has  predicted 
for  waves  originated  in  the  auroral  electrojet.  No  packets  were  registered  with  the 
period  of  oscillations  decreasing  with  time.  Tliese  results  are  consistent  with  obser¬ 
vations  of  Beggs  and  Butcher  ( 1 989). 

Relatively  small  number  of  analysed  packets  did  not  permit  to  estimate  some 
temporal  statistic  parameters  for  various  ranges  of  periods  of  oscillations.  So  we  had 
to  do  it  for  a  total  range  of  pej^ods  fi-om  T=  7  min  to  T=  50  min.  The  average  dura¬ 
tion  of  wave  packets  was  9^jpin  with  average  number  of  oscillations  of  five.  The 
maximum  duration  of  wave  p^kets  was  observed  to  be  250  min  for  a  wave  period 
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of  50  min  and  180  min  for  the, wave  period  of  20  min.  The  minimum  duration  of 
wave  packets  was  observed  to  be  50  min. 


The  frequency  -  window  procedure  was  applied  to  observations  of  midlati¬ 
tude  TIDs  at  several  heights.  This  procedure  gave  an  opportunity  to  reveal  some 
new  properties  of  TIDs.  Generally,  TIDs  are  composed  of  temporal  intervals  with 
quasi-stochastic  oscillations  of  the  electron  density  and  intervals  having  packet-like 
structures  with  quasi-monochromatic  oscillations.  Computing  cross-correlation  co¬ 
efficients  for  both  types  of  TIDs  and  applying  method  developed  at  Appendix,  we 
found  the  coherence  to  be  close  to  1  even  for  high  frequencies  TIDs  progressed  in  ■ 
the  vertical  plane  up  to  distance  80  km.  At  first  sight  it  is  not  consistent  with 
Yakovets’  et  al.  (1995)  measurements  in  which  the  fast  drop  of  the  coherence  of 
the  high  frequency  component  TIDs  was  proved.  In  fact,  there  is  not  a  contradiction. 
The  results  of  the  previous  paper  were  obtained  by  processing  the  whole  records 
which  lasted  6-8  hours  and  comprised  both  quasi-stochastic  and  packet-like  TIDs. 
Values  of  the  coherence  depended  on  a  contribution  of  wave  packets  in  the  total 
TIDs  field  —  a  ratio  of  amplitudes  of  packets  to  those  of  quasi-stochastic  fluctua¬ 
tions  and  a  ratio  of  their  lengths.  As  a  result  the  ’’average”  coherence  was  esti¬ 
mated. 

In  order  to  explain  a  difference  between  vertical  phase  velocities  of  the  wave 
packets  and  quasi-stochastic  TIDs,  one  has  to  apply  to  the  full  correlation  analysis 
(FCA)  method  developed  to  obtain  parameters  of  turbulence  from  records  registered 
by  spaced  receivers  (Briggs  et  al.  ,1950  ).  Authors  have  shown  that  the  apparent 
velocity  computed  from  a  lag  of  the  cross-correlation  function  for  records  of  spaced 
receivers  exceeded  the  true  drift  velocity  of  turbulent  structures.  A  number  of  pa¬ 
pers  were  devoted  to  this  phenomenon,  but  it  is  beyond  the  scope  of  this  paper.  Re¬ 
cent  works  of  Dewan  (1985,1991)  have  developed  the  gravity  wave  power  spectra 
model  based  on  the  turbulent-like  wave  field  picture.  This  model  seems  to  be  pres¬ 
ent  quasi  -  stochastic  TIDs  resulted  from  the  multiple-waves  superposition.  A  si¬ 
militude  of  descriptions  of  turbulence  and  the  quasi-stochastic  TIDs  field  permits  the 
FCA  method  to  be  used  for  processing  of  TIDs  records.  Although  there  is  no  full 
understanding  of  how  TIDs  parameters  can  be  related  with  results  of  an  application 
of  the  FCA  method,  a  difference  between  vertical  velocities  of  wave  packets  and 
quasi-  stochastic  TIDs  appears  to  be  explained  by  some  peculiarities  of  the  FCA 
method  mentioned  above. 

There  are  two  theoretical  models  of  a  behaviour  of  an  oscillation  period  in¬ 
side  a  wave  packet.  Francis  (1974),  considering  the  auroral  electrojet  as  a  source 
of  TIDs,  has  shown  that  after  travelling  through  the  atmosphere  up  to  F  region  the 
earth-reflected  waves  gain  the  appearance  of  a  wave  packet.  The  oscillation  fre¬ 
quency  inside  the  packet  decreases  with  time,  approaching  to  an  asymptotic  value. 
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The  asymptotic  and  average  frequencies  of  a  wave  packet  decrease  linearly  with 
growing  of  a  distance  between  a  source  of  TIDs  and  a  point  of  observations. 

An  analysis  of  the  correlation  between  magnetic  activity  and  five  occurrences 
of  packets  observed  at  Alma-Ata  with  the  growing  periods  of  oscillations  did  not 
show  evident  exceeding  of  Kp-index  for  these  days  compared  with  days  when  a  pe¬ 
riod  of  oscillations  inside  wave  packets  was  constant. 

Sources  located  in  the  troposphere  would  generate  direct  and  earth-reflected 
waves  which,  superimposing,  would  give  a  single  wave  packet  in  the  F  region.  On 
the  basis  on  the  Hines’  dispersion  relation  (1960),  Waldock  and  Jones  (1987) 
have  suggested  that  TEDs  associated  with  the  tropospheric  jet  stream  appear  in  the  F 
region  as  a  wave  packet  with  quasi-monochromatic  oscillations  whose  period  is  a 
fimction  of  the  elevation  angle  of  the  wave  ray  path  from  the  source  to  the  point  of 
observations  in  the  F  region.  Absence  of  meteorological  data  did  not  allow  us  to 
compare  them  with  TIDs  activity.  However,  as  Beggs  and  Butcher  (1989)  pointed 
out  that  even  with  data  available  it  was  not  possible  to  prove  a  relationship  between 
the  tropospheric  jet  stream  and  the  occurrence  of  TIDs  “since  TIDs  represent  an 
imperfect  tracer  of  the  passage  of  gravity  waves  in  the  thermosphere”. 

One  of  the  reason  of  phase  jumps  have  been  discussed  by  Waters  et  al. 
(1991).  They  studied  Pc  3-4  geomagnetic  pulsations  which  were  often  observed  as 
separate  wave  packets.  Numerous  records  have  shown  that  jumps  of  the  current 
phase  occur  between  close  successive  packets.  An  analysis  of  simulated  pulsations 
showed  that  phase  jumps  can  be  explained  in  terms  of  phase  variations  of  a  resul¬ 
tant  wave  formed  by  superposition  of  two  adjacent  wave  packets  having  nearby  fre¬ 
quencies  of  oscillations.  A  value  of  the  phase  jump  depends  on  the  ratio  of  ampli¬ 
tudes  of  the  packets  superimposed.  The  smaller  a  difference  of  their  amplitudes,  the 
larger  the  value  of  the  phase  jump.  If  the  envelopes  of  interfering  packets  have  equal 
amplitudes,  then  a  180®  jump  occurs  in  the  resultant  wave  packet.  At  this  moment 
the  minimum  amplitude  of  the  resultant  packet  occurs. 

Different  phase  jumps  observed  at  different  heights  at  one  series  ( Fig.  32  b  ) 
seem  to  be  caused  by  the  different  ratio  of  amplitudes  of  interfering  packets.  The 
more  value  of  the  phase  jump  took  place,  the  less  amplitude  corresponded  to  this 
moment  occurred.  Therefore  phase  jumps  observed  may  be  explained  as  a  result  of 
the  superposition  of  two  wave  packets  with  nearby  frequencies  of  oscillations.  All 
these  packets  were  quasi-monochromatic  and  hence  can  be  associated  with  the  tro¬ 
pospheric  jet  streams. 

Francis  (1974)  has  predicted  direct  and  earth-reflected  waves  from  tropo¬ 
spheric  sources.  In  order  to  explain  beating  of  these  waves,  one  have  to  suggest  the 
motkWl  of  the  jet  stream  which  would  provide  a  different  Doppler  shift  for  the  both 
w^vcs.  It  does  not  seem  to  be  possible  if  the  earth  is  plane  in  the  vicinity  of  a  region 
of  the  wave  generation  because  the  direct  and  earth-reflected  waves  would  travel 
at  the  same  angle  with  respect  to  the  horizontal  plane  to  superimpose  uvti|^q  F  re- 
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gion.  However,  a  proximity  of  Almaty  to  the  large  mountain  region  may  change  a 
condition  of  equality  of  angles  and  so  lead  to  different  Doppler  shifts. 

Different  values  of  phase  jumps  changed  phase  lags  between  adjacent 
heights,  but  it  took  place  only  in  short  intervals  of  the  order  of  a  wave  period  after 
which  the  phase  lags  became  the  same  as  in  the  previous  packets.  Tlierefore,  verti¬ 
cal  phase  velocities  have  to  be  measured  in  the  temporal  intervals  with  large  ampli¬ 
tudes  to  obtain  reliable  results.  It  is  consistent  with  criteria  used  by  JACOBSON  et  al. 
(  1995)  in  their  experiment. 

Travelling  ionospheric  disturbances  observed  over  Almaty  by  the  multifre¬ 
quency  Doppler  sounder  in  1988-1990  have  been  analysed  by  using  complex  de¬ 
modulation  method.  The  various  properties  of  two  types  of  TIDs  called  quasi¬ 
stochastic  and  wave  packets  were  studied.  A  method  of  a  calculation  of  the  cross¬ 
correlation  coefficient  for  the  wave  packets  propagating  on  the  background  of  quasi¬ 
stochastic  TIDs  was  developed.  The  high  coherence  in  the  vertical  plane  for 
’’pure”  wave  packets  was  found.  Most  wave  packets  had  quasi-monochromatic 
structure  explained  on  the  basis  of  the  wave  propagation  from  the  sources  of  TIDs 
located  in  the  troposphere  .  A  small  number  of  wave  packets  had  oscillations  with  a 
period  increased  with  time.  This  behaviour  followed  the  Francis’  (1974)  prediction 
for  earth-reflected  medium  gravity  waves  generated  in  the  auroral  region. 


APPENDIX 


Let  us  denote  by  S\{i)  oscillations  of  the  Doppler  sliift  caused  by  coherent 
fluctuations  of  the  electron  density  and  by  n{t)  oscillations  caused  by  quasi¬ 
stochastic  fluctuations  of  the  electron  density.  A  symbol  1  will  define  the  fluctua¬ 
tions  for  the  height  h\  and  a  symbol  2  will  be  used  for  the  height  h2.  Evidently  that 
quasi-stochastic  fluctuations  caused  by  numerous  sources  of  AGWs  exist  continu¬ 
ously,  so  when  coherent  packets  are  observed  the  total  values  of  the  Doppler  shift 
will  be  the  sum  of  the  coherent  and  quasi-stochastic  fluctuations: 

jc,  (t)  =  (t)  +  «,  (t)  and 

Ar,(/)  =  Sj(/)  +  n,(r) 


for  the  heights  hi  and  h2 ,  respectively. 

Mean  values  jc(/),  S(t),  and  n{t)  equal  to  0  because  of  preliminary  band-pass 


filtering  of  original  Doppler  records  and  thus  the  cross-correlation  function  may  be 
written,  as: 


[5'|(<)  +  ni(0y.S,(>  +  r)  +  «,(?  +  r)] 


(5) 
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where  overscored  bar  shows  the  mean  of  the  product  and  cr^,(7s 
are  dispersions  of  respective  values. 

Assume  that  maximum  is  centered  at  r  =  0,  and  further  we  will 

write  instead  R^^^^  (O) .  Processes  S(t)  and  «(t)  are  independent  and  so 

5'i  (/ )«2  (0  =  ‘^2  (0”i  (0  “  ®  ■  Hence 


~  - 

Jfl  X2 


(6) 


Denote  a  ratio  of  the  power  of  the  coherent  part  of  oscillations  to  the  power 
of  the  noise  as 

_2 


and  coefficient  of  correlation  as 


s,(t)^s,(0 


and 


Kn  = 

"l"2 


^  Si  ^$2 


cr« 

n,  «2 


Taking  into  account  that 


^2  2  ,  _2 


(7) 


we  obtain 


R..  =R. 


Ml 


^1  +  ql  ^ 


+  q; 


+  R 

2  '*l"2 


1 


+  q^^Ju^ 


(8) 


Here  all  values  except  are  determined  from  the  experimental  records. 

Note  that  filtered  Doppler  records  content  quasi-stochastic  temporal  parts 
which  give  an  opportunity  to  obtain  R^^ ,  and  and  temporal  parts  with 

wave  packets  which  allow  to  estimate  R^^^^ ,  and  .  Assuming  that  quasi¬ 
stochastic  Doppler  fluctuations  are  stationary  process,  one  can  think  that  values  of 
parameters  obtained  from  quasi-stCNefeastic  temporal  parts  of  records  are  valid  for 
temporal  parts  comprising  wave  pac|^^  in  order  to  compute  gI  and  R^^g^  by  using 
(7)  and  (8)  respectively. 


2.4  NONLINEAR  WAVE  -  WAVE  INTERACTIONS  OF  TRAVELLING  IONOSPHERIC 

DISTURBENCES 

Acoustic  gravity  waves  (AGWs)  propagating  in  the  neutral  atmosphere  are 
manifistated  in  the  ionosphere  as  traveling  ionospheric  disturbances  (TIDs).  It  was 
shown  (  Kaliev  et  al.,  1988)  that  TDDs  spectra  are  composed  of  a  series  of  spectral 
peaks  existing  on  the  background  of  continuous  spectra  whose  density  degreases 
toward  the  high  frequencies  edge  of  a  spectral  band.  The  line  typed  TIDs  spectra 
give  a  possibility  to  carry  out  experimental  studying  of  wave  -  wave  interactions. 

It  is  known  (Turner,  1973)  that  nonlinear  resonant  interaction  among  three 
dispersive  waves  can  occur  when  the  conditions  of  k^+k^  =  k^  and 
0,(^,)Tf  ©2(^2)  =  (o^(k^)  are  satisfied.  The  new  wave  with  frequency  co^ijk^)  originated 
into  wave  -  wave  interaction  grows  linearly.  Extensive  theoretical  investigations  of 
the  interactions  were  carry  out  for  waves  in  the  ocean  (Miller  et  al.  1986)  and  the 
atmosphere  (Dong  and  Yeh,  1988).  However,  small  experimental  studies  were  made 
for  TIDs  propagating  in  the  ionosphere.  Ma  and  Schlegel  (1993),  analyzing  data 
obtained  by  the  EISCAT  facility,  found  that  nonlinear  wave  -  wave  interactions  took 
place  for  TIDs  spectral  components  in  the  auroral  F  -  region. 

This  conclusion  was  made  on  the  basis  of  an  analysis  of  a  single  TIDs  record. 
Therefore,  to  enhance  a  statistical  reliability  of  this  conclusion  additional  studying  is 
required. 

The  purpose  of  this  section  is  to  study  nonlinear  wave  -  wave  interactions  of 
mid  -  latitude  TIDs  and  to  compare  them  with  results  obtained  by  Ma  and  Schlegel 
(1993). 

Here  we  describe  briefly  bispectral  analysis  imployed  to  estimate  a  fact  of 
nonlinear  interactions.  Suppose  x(t),  t  =  1,2  ....N,  is  a  series  of  experimental  data. 
The  series  is  subdivided  into  p  segments  each  of  length  M.  The  complex  spectrum 
of  the  j-group: 

A/  AQ\ 

f=l 

where  Wj  (t)  -  is  a  "  time  window  "  weighting  fimction  used  for  improving  of  reli¬ 
ability  of  spectral  estimates.  After  that  complex  spectra  are  summed  to  obtain  the 
spectral  average  estimate: 
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where  the  asterisk  denotes  the  complex  conjugate. 

The  segments  may  be  overlapped  or  disjoint.  If  disjoint  segments  are  used  a 
noticeable  fraction  of  the  series  will  be  lost  because  weighted  x(t)  are  close  to  zero 
near  the  edges  of  segments.  So  segments  overlapping  has  to  be  used  to  enlarge  the 
percentage  of  the  series  being  involved  into  the  S(ro)  estimates. 

The  estimate  of  the  bispectrum  is  given  as  (Nikias  and  Raghuveer,  1987) 

1  ^ 

B(,o)^,o}^)  =  -2  Xj{(o,)  •  Xj{o}^)-  +  (o^), 

P  j=\ 


The  bispectrum  gives  a  measure  of  the  multiplicative  interaction  of  the  three 
frequency  component  of  the  x(t)  spectrum.  The  normalized  bispectrum  is  the  bico¬ 
herence; 


b^{o},,o}^)  = 


P  7  =  1  P  7=1 


(11) 


If  the  component  at  /j  =  /,  +  is  generated  by  nonlinear  interaction  and  ad¬ 
ditional  unwanted  noise  is  significantly  less  of  the  interactive  components  then  the 
bicoherence  is  close  to  1.  The  necessary  condition  for  b^(co^,a)j)  to  be  close  to  I  is 
'F3=T'i+T^2,  where  'Ti  -  is  a  phase  of  the  correspondent  frequency.  We  realized  a 
program  of  bispectral  analysis  at  PC  "  Pentium".  In  order  to  check  up  its  running  an 
artificial  time  series  was  generated  as: 

x{t)  =  cos(2;^,  +  'F, )  +  cos{27rfj  +  '^2 )  "'■  cos(2;z/'3  +  Tj)  +  /;(/),  (12) 

with  frequencies  chosen  nearby  real  frequencies  of  Doppler  shift  fluctua¬ 
tions;  fi=0.025  min'\  f2=0.043  min'  and  f3=0.068  min’.  The  initial  phases  were 
satisfied  to  T^3=^Pi+'P2.  n(t)  -  is  an  additive  white  noise. 

The  series  consisted  of  600  data  points  at  1  -  min  intervals,  giving  a  total  rec¬ 
ord  length  of  10  hour.  The  Hemming  window  were  used: 

In  M  M 

W (t)  =  0.54  +  0.46  •  cos(— —  r),  r  =  -  1,0,1... 

M2  2 

To  provide  a  good  spectral  resolution  the  window  width  was  chosen  to  be 
M=120.  Segment  overlapping  was  equal  to  50%.  So,  there  were  10  segments  for  the 
estimation  of  Fourier  coefficient  according  to  (9).  Artificial  series  were  processed 
by  the  bispectral  analysis.  The  spectrum  calculated  was  composed  of  three  peaks  at 
frequetijcj^si  fi,  f2,  fs.  Fig.33  a,  b  show  the  bispectrum  and  bicoherence  of  series 
consi49rp4-  The  bispectrum  and  bicoherence  levels  are  plotted  as  counters. 


43 


The  bispectrum  and  bicoherence  reveal  the  peak  at  frequencies  fi«0.025  min  *, 
f2«0.043  min'*.  The  value  of  this  peak  for  the  bicoherence  is  close  to  1. 

Others  maxima  of  b^{a)^,(o^)  scattered  in  the  fi-ff2  domain  are  caused  by  addi¬ 
tive  noise  n(t).  The  series  (12)  is  an  idealized  model  of  interacting  TIDs  compo¬ 
nents.  Measurements  show  (section  2.3  )  that  TIDs  exist  as  packets  of  oscillations 
taken  places  during  a  fraction  of  total  records.  Of  course,  this  fact  have  to  decrease 
a  value  of  the  bicoherence. 

The  bispectral  analysis  was  applied  to  Doppler  shift  records  used  in  the  sec¬ 
tion  2.3  in  connection  with  the  study  of  packet  like  TID  structures.  Fig.  34  a,  b  show 
the  bispectrum  and  bicoherence  of  a  typical  experimental  record  obtained  in  Feb. 
23,  1989  between  10.40  and  18.40  LT  at  the  operating  frequency  of  6.7  MHz.  One 
can  see  a  chain  of  spectral  maxima  along  the  line  of  F  =  0.02  min  *.  Tables  com¬ 
prised  frequencies  participated  into  nonlinear  interactions  and  values  of  correspon¬ 
dent  bicoherences. 


Table  9 


N 

fi 

fi 

fi+fi 

1 

0.02 

0.04 

0.06 

0.6 

2 

0.02 

0.06 

0.08 

0.6 

3 

0.02 

0.11 

0.13 

0.4 

4 

0.02 

0.13 

0.15 

0.6 

5 

0.045 

0.055 

0.11 

0.6 

6 

0.08 

0.12 

0.2 

0.5 

It  is  seen  from  Table  9  that  many  interacting  frequencies  are  a  sum  of  other 
interacting  pairs.  There  seem  to  be  the  first  order  interactions  ( triplets  number  1,  5) 
and  the  high  order  interaction  (triplets  number  2,  3, 4,  6). 

Fig.  31  a,  b  show  the  bispectrum  and  bicoherence  of  an  experimental  record 
obtained  in  Feb.  12,  1989  between  10.00  and  17.00  LT  at  the  operating  frequency 
of  6.3  MHz  .  From  a  picture  it  can  be  seen  a  series  of  interacting  triplets  listed  in 
the  Table  10. 
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Table  10 


N 

fi 

fi 

fl+f2 

b'CfiA) 

1 

0.04 

0.065 

0.105 

0.6 

2 

0.035 

0.105 

0.140 

0.5 

3 

0.04 

0.135 

0.175 

0.7 

4 

0.075 

0.135  ' 

0.21 

0.7 

Similarly  the  record  obtained  in  Feb.  23,  1989  this  one  demonstrates  the  first 
and  high  order  interactions.  The  total  number  of  records  processed  was  13.  All  rec¬ 
ords  revealed  wave  -  wave  interaction  between  spectral  components  of  TIDs. 

Following  Ma  and  Schlegel  (1993)  we  can  conclude  that  existing  of  high  or¬ 
der  interaction  in  our  records  shows  an  energy  cascading  fi:om  long  periods  to 
shorter  ones.  However,  comparing  Ma  and  Schlegel  and  our  results  we  found  one 
discrepancy  deserving  to  be  considered.  Values  of  bicoherences  under  our  consid¬ 
eration  did  not  exceed  0.8  and  their  median  value  was  0.6.  At  the  same  time  Ma 
and  Schlegel  (1993)  estimated  bicoherences  to  be  0.996,  0.98  and  0.98  for  different 
interacting  pairs.  Such  high  values  of  bicoherences  for  experimental  records  are  un¬ 
usual  and  could  be  only  expected  for  artificially  generated  data  with  a  low  level  of 
the  additional  noise,  where  the  conditions  (coi+(D2=g)2,  M^i+V2=M^3)  of  a  resonant  in¬ 
teraction  between  the  full  time  existing  fi’equencies  are  valid. 

These  high  values  of  bicoherences  seem  to  be  caused  by  a  relatively  short 
length  (5  hours)  of  the  record  processed  by  Ma  and  Schlegel  (1993).  In  spite  of  dis¬ 
crepancy  considered  above  following  conclusions  seem  to  be  evident: 

1 .  Nonlinear  wave  -  wave  interactions  are  common  feature  of  TIDs  propaga¬ 
tion  in  the  F  -  region  of  the  ionosphere. 

2.  Existing  of  the  high  order  nonlinear  interactions  appears  to  prove  an  energy 
cascading  firom  long  periods  to  shorter  ones. 
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Fig. 18.  Plots  of  virtual  height  variations  in  the  night  25-26.12.1997 
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Fig. 21.  Plots  of  virtual  height  variations  in  the  night  5-4.04,1998. 
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Fig.22.  Distributions  of  wave  periods,  vertical  phase  speeds,  wavelength  and  relative 
amplitudes  for  waves  analyzed. 
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Fig.  23.  Relative  amplitude  versus  Kp-index  (a),  averaged  relative  amplitude  versus 
month  of  observations  (b). 
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Pig. 24.  Detrended  variatios  of  normalized  virtual  height 
at  the  frequency  f^  (a),  the  frequency  f^  (b)  and  their 
normalized  power  for  Nov. 27,1997. 
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operating  frequencies  f^=5.4  MHz,  f2+5.7  MHz,  fj=6,3  MHz,  f^=6,7  MHz,  f^= 
7.4  Mz  and  passed  through  the, filter  with  the  central  frequency  F=1,5  cph 
The  vertical  her  denotes  the  0,5  Hz  scale  of  the  Dqppler  shift  variations. 


shift  variations  nlotted  in  j'ia*.  24 
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shift  variations  plotted  in  i'ig,27 
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Fig, 33.  Bispectrum  (a)  and  bicoherence  (b)  for 
artificial  series. 
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Fig .3^  «Bispectrum  (a)  and  bicog_erence  (b)  for  the  experimental 
record  obtained  in  Feb. 23, 1989* 
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ITig,  55.  Bispectrum  (a)  and  Mcogerence  (b)  for  the  experimental 
record  obtained  in  Feb. 12, 1989. 
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lEI.  Simulation  of  Planetary  Waves  and  Their  Influence  on  the  Zonally 
Averaged  Circulation  in  the  Middle  Atmosphere 

Paper  was  presented  by  A.I.Pogoreltsev  as  an  invited  report  at  the  International  Symposium 
on  Dynamics  and  Structure  of  the  Mesopause  Region  (March  16-20,  1998,  Kyoto,  Japan)  and 
will  be  submitted  for  publication  in  Special  Issue  of  Earth,  Planets  and  Space  in  June  1998. 

3.1.  Introduction. 

Analysis  of  the  results  obtained  by  the  Wind  Imaging  Interferometer  (WINDII)  and  High  Res¬ 
olution  Doppler  Imager  (HRDI)  on  the  Upper  Atmosphere  Research  Satellite  (UARS)  shows  a 
presence  of  the  strong  large-scale  variations  of  the  emission  rates  and  winds  in  the  mesosphere 
and  lower  thermosphere  (Shepherd  et  al.,  1993;  Smith,  1995),  which  were  interpreted  as  man¬ 
ifestation  of  the  planetary  waves  (PWs).  Besides,  the  main  purpose  of  the  PSMOS  project 
is  the  establishment  of  a  ground-based  mesopause  observing  system  for  the  investigation  of 
planetary  scale  disturbances,  and  the  purpose  of  present  report  is  to  show  some  results  of  the 
simulation  of  PWs  in  the  middle  atmosphere  and  to  discuss  a  possibility  of  the  observation 
of  these  waves  by  satellite  and  ground-based  equipment.  The  comparison  of  numerical  results 
with  the  experimental  data  will  be  made  and  the  contribution  of  PWs  to  the  zonally  averaged 
dynamics  of  mesopause  region  will  be  considered. 

3.2.  Model  Formulation. 

The  initial  set  of  the  equations  describing  the  disturbances  and  zonally  averaged  state  can 


be  written  as  follows  (Boyd,  1976;  Andrews,  1985) 

{u  -  c)u\l (a  cos  (p)  -  /iv'  +  u^w'  -f-  </>^/(o  cos  ip)  =  X\  (1) 

(u  -  c)v'J {a  cos  ip)  +  /2u'  +  ^'^/ a  =  Y\  (2) 

if>',  =  Re'exp{-KzlH)/H,  (3) 

u'J (a  cos  ip)  +  (cos  <pv')^/ {a  cos  ip)  -f  {pow')^lpQ  =  0,  (4) 

(u  -  c)^5^/(a  cos  (/?)  -H  O^v'l a  +  =  Q',  (5) 

-  fiv*  -t-  UzW*  =  X  +  V‘F/{poa  cos  ip),  (6) 

{f  +  utanip/a)u  +  —  G,  (7) 

=  RO  Gxp{—Kzf  H)/H,  .  (8) 

(cos  ipv*)^/ {a  cos  ip)  4-  {pow*)Jpo  =  0,  (9) 

Ot  +  e^vya  -h  e,w*  =  q  -  {poxe-  Vd'  ih)zlpo.  (lo) 
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Here  A  is  longitude,  (p  is  latitude,  z  = -H  ln(p/ps)  is  the  vertical  coordinate,  where p  is  pressure, 
Pg  =  const  is  a  standard  reference  pressure,  H  =  const  is  the  scale  height;  u,  v,  and  w  are 
the  corresponding  velocities,  0  denotes  the  geopotential,  and  6  the  potential  temperature; 
R  is  the  gas  constant  for  dry  air,  k  =  Rjop,  Cp  is  the  specific  heat  at  constant  pressure; 
X,  Y  are  the  unspecified  horizontal  components  of  friction  or  other  nonconservative  mechanical 
forcing,  Q  is  the  diabatic  heating  term;  po{z)  =  psexp{-zfH)  is  a  reference  density  profile. 
Pa  =  const,  /i  =  /  +  utancp/o  —  u^ja,  /2  =  /  +  2utan(p/a,  /  =  20  sin  9?  is  the  Coriolis 
parameter;  0  and  a  are  the  Earth’s  rotation  rate  and  radius;  F  =  (0,  is  the  EP  flux 

vector,  the  components  of  which  being  given  by  the  following  relations 


=  Pqo  cos  ip{uz  v'6'  I6z—  u'v'),  (11) 


=  pacos(p{fi  v'9'  I0z~  u'w),  (12) 

G  represents  all  the  terms  that  lead  to  departures  from  the  gradient  wind  balance  between 
and  u;  v*  and  w*  define  the  residual  meridional  circulation  velocities,  which  are  related  with 
Eulerian-mean  ones  as  follows 

V*  =  V  —  {po  v'9'  I9z)zIpo,  (13) 

w*  =  w  +  (cos  ip  v'9'  cos  (/?).  (14) 

In  equations  (I)-(IO)  an  overbar  denotes  the  zonally  averaged  values,  the  disturbance  is  denoted 
by  a  prime;  the  subscripts  A,  tp,  z,  t  denote  corresponding  derivatives. 

We  assume  the  disturbances  in  equations  (I)-(IO)  to  be  a  sum  of  zonal  harmonics  ,  i.e., 
q'{X,<p,z,t)  =  Fei:m9m(v’>^)exp[f("iA  -wt)],  where  q'^  is  the  complex  amplitude  and  m  is 
the  zonal  wavenumber.  It  is  necessary  to  note  that  equations  (l)-(5)  are  the  linearized  set  of 
equations  for  disturbances,  and  that  they  apply  for  a  single  zonal  harmonic  with  phase  speed  c. 
The  formulation  of  the  boundary- value  problem  for  simulation  the  latitude  height  structure  of 
the  planetary  scale  disturbances  and  the  numerical  method  of  the  solution  were  described  in  the 
paper  Pogoreltsev  and  Sukhanova  (1993),  and  here  we  will  only  present  and  discuss  the  results 
of  the  simulations  with  this  model.  However,  it  is  necessary  to  note  here  that  the  simplification 
of  the  initial  set  of  equations  for  disturbances  based  on  the  scale  analysis  [by  assuming  a  large 
Richardson  and  a  small  Rossby  numbers  (Pogoreltsev  and  Sukhanova,  1993)]  are  not  justified 
in  the  case  of  inertia-gravity  PWs  [for  instance,  2-  and  5-day  travelling  PW],  and  in  present 
paper  we  will  consider  the  solution  of  exact  primitive  set  of  equations  (10-(5). 

3.3.  The  Background  State. 

The  numerical  simulations  of  the  global  structure  of  planetary  waves  were  carried  out  for 
the  climatic  conditions  of  the  Northern  Hemisphere  winter.  The  background  temperature  was 
taken  from  the  MSIS-90  model  (Hedin,  1991).  The  vertical  profile  of  Tq{z)  is  shown  in  Fig.  1. 
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The  simulations  that  we  discuss  were  performed  for  the  atmosphere  with  scale  height  H—7  km 
using  the  vertical  step  Az  =  H/A. 

We  used  two  models  of  the  background  wind  for  January  conditions 

I.  HWM-93  model  (Hedin  et  ah,  1996); 

II.  In  altitude  region  between  90  and  120  km  we  used  an  empirical  model  of  zonal  wind  based 
on  the  WINDII  data  (Wang  et  ah,  1997). 

The  background  wind  distributions  in  the  HWM-93  and  WINDII  models  are  shown  in  Fig. 

2.3.  One  can  see  that  in  these  models  the  westerlies  over  equator  in  the  mesosphere  are  absent. 
However,  the  direct  measurements  by  the  HRDI  for  January  1992  through  June  1993  show  the 
westerlies  over  equator  up  to  the  90  km  in  winter  (Lieberman  et  ah,  1993).  To  investigate  the 
effect  of  the  mesospheric  westerlies  over  equator  on  the  propagation  of  the  PWs,  we  used  the 
MODIFIED  WINDII  model,  which  was  defined  as  follows 

u  =  «(WINDII)  30exp[-(^^)'*]  exp[-(^)2].  (15) 

The  distribution  of  the  background  wind  for  the  MODIFIED  WINDII  model  is  shown  in  Fig.  4. 

3.4.  Dissipation  Parameters. 

The  dissipative  processes  were  parameterized  in  the  form  of  Rayleigh  friction  and  Newtonian 
cooling  with  coefficients  I3r{z,  (f)  a.nd  aN{z).  To  avoid  the  numerical  instability  of  the  solution 
near  the  critical  surfaces,  we  introduce  a  localized  effective  Rayleigh  friction  coefficient  in  these 
regions 

(3r{z,  v?)  =  (3ro  +  {mpRm  -  Pro)  exp[-( — - — ^ — )^],  (16) 

a  cos  ^fjRrn 

where  Pro  =  10“®s~^,  PRm  =  3  •  10“®s“h  For  Newtonian  cooling  we  used  the  parameterization 
of  damping  rate  in  the  cool-to-space  limit  given  by  Zhu  (1993).  Above  120  km  the  infrared 
radiation  resulting  from  5.3  uva.  nitric  oxide  band  plays  the  main  role  in  the  lower  thermosphere 
heat  balance,  and  we  used  the  formula  given  by  Koskarts  (1980)  to  calculate  the  damping  rate 
in  the  thermosphere.  The  resulting  Newtonian  cooling  coefficient  is  presented  in  Fig.  5. 

3.5.  Boundary  Conditions. 

The  amplitudes  and  phases  of  geopotential  height  perturbations  in  stationary  PW  with  m=l 
at  the  level  18.55  mb  {z=29>  km)  taken  for  January  from  CIRA-86  climatic  model  (Barnett  and 
Corney,  1985)  were  set  as  a  low  boundary  condition.  The  perturbation  of  the  geopotential  at 
the  upper  boundary  which  is  set  at  140  km  was  taken  to  be  equal  zero. 


Background  wind  (WINDII,  January) 


Figure  3:  As  in  Fig.  2,  but  for  WINDII  model 
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3.6.  The  Influence  of  Background  Wind  Distribution  on  SPW. 

The  results  of  simulation  of  the  stationary  planetary  wave  (SPW)  with  zonal  wave  number 
m=l  in  a  form  of  isolines  of  amplitude  of  the  geopotential  height  perturbation  are  shown  in 
Fig.  6-8  for  HWM-93,  WINDII,  and  MODIFIED  models  of  the  background  wind.  One  can 
see  that  in  the  cases  of  the  HWM-93  and  WINDII  models  SPW  propagates  mainly  into  the 
high  latitude  region  of  the  Northern  Hemisphere  and  then  dissipates  at  the  heights  of  the  lower 
thermosphere  (see  Fig.  6,7).  In  the  case  of  MODIFIED  model  there  is  a  penetration  of  SPW 
into  the  Southern  Hemisphere  across  the  equator  through  the  mesospheric  westerlies.  In  result 
the  amplitudes  of  the  geopotential  height  perturbations  in  the  SPW  with  m=l  are  comparable 
in  the  lower  thermosphere  of  the  winter  and  summer  Hemispheres. 

3.-7.  Comparison  of  Computations  and  Observations. 

The  vertical  structure  of  the  geopotential  height  perturbation  in  SPW  with  m=l  at  the  lat¬ 
itude  60*^N  is  shown  in  Fig.  9.  Dotted  lines  show  the  vertical  profiles  of  the  amplitude  and 
phase  taken  from  empirical  model  of  SPW  (Barnett  and  Cornney,  1985).  The  comparison  of 
the  results  of  simulation  (solid  lines)  with  observations  shows  good  agreement. 

The  latitudinal  dependencies  of  the  amplitude  and  phase  of  the  meridional  wind  perturbation 
in  SPW  with  zonal  wave  number  m=l  at  the  altitude  z=90  km  are  shown  in  Fig.  10.  Dotted 
lines  show  the  results  obtained  from  WINDII  data  averaged  over  period  1992-1996  (Wang  D.Y., 
private  communication).  Solid  lines  show  the  results  of  simulation  with  the  MODIFIED  back¬ 
ground  wind  model.  One  can  see  that  the  latitudinal  structure  of  the  amplitude  does  not 
show  good  agreement  between  \VINDII  observations  and  solution.  Perhaps  the  reason  for  this 
discrepancy  is  the  strong  latitudinal  smoothing  of  the  background  wind  due  to  including  only 
limited  set  of  the  vector  spherical  harmonics  in  the  analytical  empirical  horizontal  wind  models 
(Hedin  et  ah,  1991;  Wang  et  ah,  1997).  However,  the  latitudinal  structure  of  the  phase  show 
good  agreement  between  WINDII  observations  and  solution.  It  is  necessary  to  note  that  west¬ 
ward  slope  of  the  phase,  when  we  move  from  the  Northern  to  Southern  Hemisphere,  indicates 
the  cross  equatorial  SPW  propagation,  and  this  result  supports  our  assumption  that  in  winter 
there  are  the  westerlies  over  the  equator  in  the  mesosphere. 


-90  -60  -30  0  30  60  90 

Latitude  (deg) 


Figure  6:  Latitude-height  section  of  amplitude  (dam)  for  the  perturbation  of  the  geopotential 
height  in  SPW  with  m=l  and  background  wind  from  HWM-93  model. 
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Figure  9:  The  vertical  profiles  of  the  amplitude  and  phase  (the  position  of  ridges)  of  the 
geopotential  height  perturbation  in  SPW  with  m=l  at  the  latitude  60°N.  Dotted  lines  show 
the  results  from  empirical  model  CIRA-85. 
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Meridional  wind,  m=l  (z=90  km) 


Latitude  (deg)  Latitude  (deg) 


Figure  10:  Latitudinal  dependencies  of  the  amplitude  and  phase  of  the  meridional  wind  per¬ 
turbation  in  SPW  with  m=l  at  the  altitude  z=90  km.  Dotted  lines  show  the  results  obtained 
from  WINDII  data  averaged  over  period  1992-1996  (Wang  D.Y.,  private  communication). 
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3.8.  Responses  of  the  Atmosphere  to  Lower  Boundary  Forcing. 

The  model  of  the  global  structure  of  SPW  (Pogorltsev  and  Sukhanova,  1993)  was  modified 
to  simulate  the  spectral  atmospheric  response  to  the  symmetric  and  asymmetric  forcing  by 
the  perturbation  of  the  vertical  velocity  at  the  lower  boundary.  This  problem  is  similar  to  the 
problem  which  was  considered  by  Schoeberl  and  Clark  (1980).  However,  there  are  some  reasons 
to  reconsider  this  question  once  more.  Schoeberl  and  Clark  (1980)  used  the  quasi-geostrophic 
vorticity  equation  for  the  simulation  of  the  structure  of  PW.  Kirushov  (1988)  investigated 
the  eigenvalues  and  eigenfunctions  of  the  latitudinal  structure  operator  of  PW  under  quasi- 
geostrophic  approximation,  and  it  has  been  shown  that  the  spectrum  becomes  quasi-confluent 
and  symmetric  modes  coincide  with  asymmetric  ones  with  the  exception  of  the  sign  in  an¬ 
other  Hemisphere  when  the  Rayleigh  friction  coefficient  tends  to  zero.  In  other  words,  we  have 
the  .complete  set  of  the  eigenfunctions  in  each  Hemisphere  separately,  and  the  Northern  and 
Southern  Hemispheres  do  not  interact  in  this  case.  In  result  low  frequency  responses  of  the  at¬ 
mosphere  to  the  symmetric  and  asymmetric  forcing  are  very  similar  (see  Fig.  2  in  Schoeberl  and 
Clark,  1980).  In  the  case  of  primitive  set  of  equations  the  responses  become  noticeably  different. 
This  problem  was  investigated  by  Salby  (1981),  and  it  has  been  shown  that  the  background 
wind  shears  cause  the  spectral  broadening  and  shifting  of  the  resonant  responses.  However,  it 
seems  that  strong  spectral  broadening  obtained  by  Salby  (1981)  is  the  consequence  of  the  low 
latitudinal  resolution  and  small  dissipation.  Lindzen  and  Blake  (1972)  showed  that  realistic 
dissipation  mechanisms  suppress  the  atmospheric  response,  and  that  the  effect  increases  with 
the  period  of  the  wave. 

To  investigate  the  sensitivity  of  the  atmospheric  response  to  lower  boundary  forcing  we  specified 
the  following  boundary  conditions.  For  the  symmetric  forcing  we  used  w*{z  =  0)  =  wScos2i?, 
and  for  asymmetric  forcing  we  used  w*{z  =  0)  =  iDqCosi?,  where  w*  =  w*2Qaf gH  cos"'  ip, 
w*  =  w'  -  is  the  analog  of  the  vertical  velocity  in  the  geometric  coordinates,  t?  =  7r/2  -  (p 
is  co-latitude,  and  wq  =  10"®.  It  is  necessary  to  note  here  that  all  calculations  connected  with 
traveling  PWs  were  made  with  the  HWM-93  model  of  the  background  wind  for  January  con¬ 
ditions. 

The  amplitude  response  of  the  model  to  the  lower  boundary  forcing  for  westward  PWs  with 
m=l  as  function  of  the  period  are  shown  in  Fig.  11.  Here  we  define  the  response  as  follows 

Response  = —^  1  exp(-2,72ff),  (17) 

where  is  a  number  of  grid  points.  One  can  see  from  this  figure  that  there  are  well  defined 
resonant  responses  at  the  periods  about  5  and  16  days  for  symmetric,  and  at  the  period  about 
10  days  for  asymmetric  forcing.  The  latitudinal  structures  of  the  geopotential  height  pertur¬ 
bations  in  the  10  and  16  days  PWs  at  the  altitude  z=7  km  are  shown  in  Fig.  12.  This  figure 
shows  that  these  waves  can  be  identified  with  first  asymmetric  (1,2)  and  second  symmetric 


91 


Geopotential  height,  10-day  PW,  z=7  km 
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Figure  12:  The  latitudinal  structure  of  the  geopotential  height  perturbations  in  the  10-  and 
16-days  PWs  at  altitude  7  km. 
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(1,3)  normal  atmospheric  modes  in  the  nomenclature  of  Longuet-Higgins  (1968).  By  analogy 
the  5-days  PW  can  be  identified  with  first  symmetric  mode  (1,1).  The  amplitude  response  of 
the  model  to  the  lower  boundary  forcing  for  the  westward  PW  with  m=2  is  shown  in  Fig.  13. 
It  is  easy  to  identify  the  maxima  in  this  figure  with  first  symmetric  (2,1)  and  first  asymmetric 
(2,2)  normal  modes  with  the  periods  about  4  and  7  days,  respectively.  The  second  symmetric 
mode  for  PW  with  m=2  does  not  develop  due  to  strong  dissipation. 

To  reproduce  the  realistic  conditions,  we  calibrate  the  solution  in  our  model  using  the  mea¬ 
sured  geopotential  height  amplitude  of  16-days  PW  in  stratosphere  averaged  over  January 
10-22,  1979  -  the  maximum  is  about  100  dam  (Smith,  1985).  The  latitude-altitude  structure 
of  the  geopotential  height  amplitude  for  16-days  PW  is  shown  in  Fig.  14.  One  can  see  that 
this  wave  propagates  from  winter  stratosphere  into  the  summer  lower  thermosphere  across  the 
equator.  We  calibrate  10  days  PW  solution  by  similar  way,  assuming  that  in  the  mesopause 
region  this  wave  has  the  amplitude  comparable  with  the  amplitude  of  16-days  PW,  and  the 
latutude-altitude  structure  of  10-days  PW  is  shown  in  Fig.  15.  The  latitudinal  structures  of  the 
geopotential  height  perturbations  in  10  and  16  days  PWs  at  the  altitude  z=96  km  are  shown  in 
Fig.  16.  One  can  see  from  this  figure  that  the  latitudinal  changes  in  these  waves  are  very  sim¬ 
ilar,  and  there  are  westward  slopes  of  phases  from  Northern  to  Southern  Hemisphere  for  both 
of  waves.  These  slopes  indicate  the  propagation  of  PW  from  winter  into  summer  Hemisphere 
across  the  equator.  We  cannot  identify  10  and  16  days  PWs  in  the  mesopause  region  with  the 
normal  atmospheric  modes  as,  for  instance,  in  the  troposphere  (see  Fig.  12),  and  it  is  possible 
to  make  a  conclusion  that  the  latitudinal  structure  of  these  waves  in  the  middle  and  upper 
atmosphere  depends  mainly  on  the  distribution  of  the  background  wind,  and  does  not  depend 
on  the  latitudinal  structure  of  the  lower  boundary  source.  The  structures  of  the  amplitude  of 
zonal  wind,  meridional  wind,  and  temperature  perturbations  due  to  16  days  PW  are  shown  in 
Fig.  17-19.  These  figures  show  that  there  exist  substantial  (drl5  m/s)  variations  of  the  zonal 
wind  in  the  mesopause  region  conditioned  by  16  days  wave.  The  variations  of  the  meridional 
wind  are  smaller  by  a  factor  of  2.  The  variations  of  the  temperature  at  the  middle  latitudes  of 
both  Hemispheres  can  reach  ±  5  K. 

The  response  of  the  model  to  asymmetric  forcing  for  westward  PWs  with  m=3  is  shown  in 
Fig.  20.  The  maximum  in  this  figure  is  situated  at  the  period  52.5  hours,  and  this  PW  can  be 
identified  with  well  known  third  Rossby-gravity  normal  atmospheric  mode  (3,0),  the  so-called 
2-days  wave.  To  reproduce  the  realistic  conditions,  we  calibrate  the  solution  using  the  observed 
by  WINDII  during  January  1993  magnitudes  of  the  horizontal  wind  variations  due  to  2  days 
wave  (Ward  et  ah,  1996).  The  latitudinal  structures  of  the  amplitude  of  zonal  and  meridional 
wind  perturbations  in  2  days  PW  at  the  altitude  z=  96  km  are  shown  in  Fig.  21,22.  Dashed 
lines  show  the  results  obtained  from  the  analysis  of  WINDII  data,  solid  lines  show  the  numerical 
solution.  One  can  see  that  latitudinal  structure  of  the  amplitude  of  the  zonal  wind  does  not 
show  good  agreement  between  WINDII  results  and  solution.  The  latitudinal  structure  of  the 
PW,  and  especially  the  perturbation  of  the^onal  wind,  is  strongly  sensitive  to  the  latitudinal 
changes  in  the  background  wind,  and  one  of  the  possible  reasons  of  the  discrepancy  between 
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Figure  16:  The  latitudinal  structure  of  the  geopotential  height  perturbations  in  the  10-  and 
16-days  PWs  at  altitude  96  km. 
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Figure  20:  Response  to  the  asymmetric  forcing  for  the  westward  PWs  with  m=3. 
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Zonal  wind 

2— day  PW,  z=96  km 


Figure  21:  Latitudinal  dependencies  of  the  amplitude  and  phase  of  the  zonal  wind  perturbation 
in  2-days  PW  with  zonal  wave  number  m=3  at  the  altitude  ^=96  km.  Dashed  lines  show  the 
results  obtained  from  analysis  of  WINDII  data 


Meridional  wind 
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Figure  22:  As  in  Fig.  21,  but  for  the  perturbation  of  meridional  wind. 
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calculated  and  observed  latitudinal  variations  of  the  zonal  wind  is  the  smoothing  of  the  back¬ 
ground  wind  in  the  HWM-93  model.  The  latitudinal  structure  of  the  meridional  wind  shows 
good  agreement  between  the  WINDII  data  and  solution.  The  latitude-altitude  structure  of  the 
geopotential  height  amplitude  in  2  days  PW  is  shown  in  Fig.  23.  One  can  see  that  2  days 
wave  propagates  from  the  equatorial  stratosphere  into  the  lower  thermosphere  of  the  Southern 
(summer)  Hemisphere.  The  calculated  structures  of  the  amplitude  of  zonal  wind,  meridional 
wind,  and  temperature  perturbations  due  to  2  days  PW  are  shown  in  Fig.  24-26. 

3.9.  The  Wave-Mean  Flow  Interaction. 

The  generalized  Eliassen-Palm  theorem  states  that  F  [see  (6)]  appears  as  the  flux  in  the  wave- 
activity  density  conservation  law,  which  can  be  written  for  waves  of  small  amplitude  a  as  follows 
(Andrews  and  McIntyre,  1976,1978) 

dA/dt  +  V-F  =  D  +  0{a%  (18) 

where  A,  D  and  F  are  mean  O(a^)  functions  of  wave  quantities  and  D  vanishes  for  conservative 
(i.e.  adiabatic,  frictionless)  waves.  Note  that  in  the  case  of  the  steady  waves  V-F  =  D.  Thus, 
for  the  steady  non-conservative  waves  D  ^  0  and  defines  the  acceleration  of  the  zonally  mean 
flow  by  the  waves. 

Andrews  (1987)  obtained  the  simple  explicit  expressions  for  the  wave-activity  density  A  and 
non-conservative  source  of  wave-activity  D,  which  involve  Ertel’s  potential  vorticity,  and  include 
only  Eulerian  quantities 

A  =  poa  cos  <p{Ai  -I-  A2  -f-  A3),  (19) 

Ai  =  ^Po¥^lJ,  A2  =  \po{uJpo)z'^  fOl  = 

D  =  Poa  cos  (p{Di  +  D2  +  D3),  (20) 


Di  =  po  H'S'  / J,  D2  =  poiujpo).  9'Q'  D3  =  -K  u',Q')/e,, 

Ii'  =  p'-d'pje,,  ^'  =  R'-Q'PJ9,,  J=-^^  =  -M-pJ^), 

d  Zj  0t 

Ertel’s  potential  vorticity  P  is  given  by 

=  +  +  (21) 

'■  acos<^  acos(/?J  acos(p  a 


and  Ertel’s  theorem  states  that 


DP/Dt  =  R, 


(22) 
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Figure  23:  Latitude-height  section  of  amplitude  (dam)  for  the  perturbation  of  the  geopotential 
height  in  2-days  PW. 


Meridional  wind  (m/s),  2— day  PW 


Figure  25:  As  in  Fig.  23,  but  for  the  perturbation  of  meridional  wind 
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where  DjDt  is  the  material  derivative,  and  R  represents  non-conservative  effects 
PqR  =  {acos^p)~^\-{YzO-\rv^Q)x^\■[{X^9^-u^Q)coslp\^^{■ 

-h{(acos(/5)“^[-(Xcos(^),^-l-yA]^+ ^  H - ^ — 1(5}  •  (23) 

t'  7"/  I  V  T  j'f'  j  a  cos  if  acosyjJ  ■'* 

The  expressions  for  P'  and  R'  are  written  as  follows 

/>oP'  =  -^[-(«'  cos  + t-g (24) 
a  cos  ip  a  a 

PqB!  =  — — [  —  {X'  cos  -f-  V^]  -f  fiQ'z  -f  - — -  H-  ^ (25) 
a  cos  p  a  a 

The  divergence  of  Eliassen-Palm  flux  can  also  be  calculated  using  the  finite  differences.  Note 
that  in  {p^  z)  coordinates 

VF  =  cos v?)J(a cosy?) (26) 


The  latitude-height  section  of  the  the  EP  flux  divergence  [the  force  per  unit  mass  on  the  mean 
flow  due  to  .PWs  dissipation]  for  SPW  with  m=l  are  shown  in  Fig.  27.  Our  calculations  show 
that  in  the  case  of  PWs  |Ai|  >  \A2\,  [Asl  (jFij  >  ID2I,  [Fs]),  and  as  a  consequence  of  this  fact, 
the  sign  of  the  wave-activity  density  (and  EP  divergence)  is  determined  by  the  sign  of  J.  The 
sign  of  J  is  mainly  determined  by  the  sign  of  the  latitudinal  gradient  of  the  background  Ertel’s 
potential  vorticity.  As  a  result,  in  the  region  of  the  negative  latitudinal  gradient  of  the  Ertel’s 
potential  vorticity  A  <  0  and  F  >  0,  i.e.,  we  have  the  eastward  acceleration  of  the  zonally 
averaged  flow. 


The  magnitudes  of  the  accelerations  due  to  the  traveling  PWs  are  comparable  with  the  ones 
due  to  SPW  with  m=l.  However,  if  we  add  together  the  effects  of  different  westward  traveling 
PWs,  the  accelerations  due  to  these  waves  become  more  significant.  In  Fig.  28  the  force  per 
unit  mass  on  the  mean  flow  due  to  2-,5-,10-,  and  16-days  waves  are  shown.  One  can  see  from 
this  figure  that  in  the  mesopause  region  these  accelerations  can  reach  the  magnitude  about  50 
ms~^/day.  This  value  is  comparable  with  the  accelerations  due  to  the  internal  gravity  waves 
and  atmospheric  tides  breakdown  and/or  saturation  (Lindzen,  1981). 


Altitude  (kra) 


EP  flux  divergence,  SPW,  m=l 
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Figure  27:  V-F/(poacos(/?)  -  the  force  per  unit  mass  on  the  mean  flow  due  to  SPW  (ni=l) 
dissipation.  Contour  interval  is  2  m/s,  solid  contours  are  positive,  dashed  negative  and  zero 
contour  is  removed. 


EP  flux  divergence,  2-.5-,10-,  and  16-day  PW 
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Figure  28:  S7-F/{poa  cos  cp)  -  the  force  per  unit  mass  on  the  mean  flow  due  to  westward  traveling 
2-, 5-, 10-,16-days  waves.  Contour  interval  is  4  m/s,  solid  contours  are  positive,  dashed  negative 
and  zero  contour  is  removed. 


3.10.  Main  Results 


The  results  of  simulations  with  the  numerical  model  of  the  SPW  and  traveling  PWs  struc¬ 
ture  can  be  formulated  as  follows 

•  There  is  a  good  qualitative  agreement  between  the  results  of  simulation  of  the  planetary 
waves  and  the  available  experimental  data. 

•  The  latitudinal  structures  of  PWs  in  the  middle  atmosphere  are  strongly  dependent  on 
the  distribution  of  the  background  wind,  and  to  provide  the  detailed  comparison  of  the 
results  of  simulation  with  the  characteristics  of  PWs  obtained  from  experimental  data,  it 
is  necessary  to  use  (and  develop)  the  background  wind  models  taking  into  consideration 
the  local  latitudinal  variations  of  the  zonally  averaged  wind. 

f 

•  There  is  the  experimental  support  (westward  slope  of  the  phase)  of  the  cross  equatorial 
propagation  pf  the  stationary  PWs  from  the  winter  stratosphere  into  the  summer  lower 
thermosphere  through  the  westerlies  in  the  low  latitude  mesosphere,  and  the  question 
about  the  existance  of  eastward  zonal  wind  in  this  region  is  crucial  for  our  understanding 
of  the  dinamical  interaction  between  Hemispheres. 

•  Planetary  waves  can  provide  the  substancial  accelerations  of  the  zonally  averaged  flow 
in  the  mesosphere  and  lower  thermosphere  so  that  to  simulate  the  general  circulation  of 
these  regions  it  is  necessary  to  develop  three  dimensional  models. 

•  The  magnitude  of  the  perturbations  of  the  meteorological  parameters  due  to  presence  of 
PWs  in  the  mesopause  region  is  sufficient  to  be  registrated  by  satellite  and  ground-based 
instruments. 
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